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SYNOPSIS 


The results of an investigation undertaken to record and 

interpret the absorption and fluorescence spectra of Eu^ + doped 

in LiYF 4 single crystal are presented in this thesis. The 

host crystal LiYF^ belongs to the tetragonal scheelite family. 

3+ 3+ 

The Eu ion substitutes for Y at a site of S^ point symmetry. 

A sample of LiYF 4 :Eu^ + , of dimensions 3.8 mm x 3.1 mm x 
2.1 mm with c-axis alono the 3.1 mm edge wa s obtained from 
Dr L.E. Erickson (NRC, Ottawa, Canada). A glass dewar employ- 
ing cold finger was fabricated for the low temperature spec- 
troscopic measurements. The absorption spectrum was obtained 
using a water cooled 1000W tungsten-halogen lamp as a white 
light source and the fluorescence spectra were recorded using 
a Coherent Ar + Laser and a Coherent ring dye laser (with R 6G 
dye). The signals transmitted or emitted by the sample were 
detected by a PMT attached to a double grating monochromator 
and then amplified by a lock-in-amplifier. These amplified 

signals were recorded on an x-t recorder. The excitation 

5 7 

spectrum was obtained by monitoring the j D q (Q)- F^CT^) 
fluorescence and scanning the dye laser within the spectral 
range of R 6G dye. The fluorescence and absorption spectra 
were obtained in %, c and axial modes at liquid nitrogen tem- 
perature. None of the Ar + laser lines is in resonance with 
any of the Eu^ + absorption transitions, but LiYF 4 :Eu^ + 


exhibits 



xi 


fluorescence when excited with any of the Ar + laser lines. 

The 457.9 nm and 514.5 nm lines produce fairly strong 

fluorescence in this sample. The selective excitation of 

5 

the D q level was obtained using the dye laser lasing at 
16938 cm ^ (which corresponds to the wavenumber interval for 
the 7 F^ ( r 3 ) “ transition) and at 17630 cm""'*' (excita- 

tion through a phonon assisted process). The excitation 

5 

spectrum was helpful in ascertaining the position of -the D q 
level and also in identifying some of the interfering Raman 
and vibronic lines. Initial assignment to the crystal field 
levels within a given multiplet was made on the basis of 
symmetry considerations and the observed polarizations in the 
absorption and fluorescence spectra. This assignment was 
later found to be consistent with our calculations. All the 
crystal field energy levels belonging to the ^Fj(J=0,l ,2,3,4) 
and “*Dj (J=0,l,2) multiplets were identified, while the 
structure of the 7 F^ and multiplets could be established 

only partially. Experimental free ion level positions were 
obtained by taking the centres of gravity of the crystal field 
levels belonging to a given multiplet. 

A detailed computational analysis of the data has been 
attempted. The free ion energy matrices were constructed 
and diagonalized to obtain the free ion energy levels and 
the wavefunc tions in the intermediate coupling scheme. The 
angular parts of the matrix elements were calculated exactly 



xii 


and the radial quantities were treated as adjustable para- 
meters. The crystal field calculations were carried out with 
and without J-mixing effects using intermediate coupling 
wavef unc tions as zero-order wavefunc tions . Here again radial 
parts v/ere parametrized in an attempt to match the experi- 
mental and theoretical results. The best possible set of 
parameters obtained within the computational constraints are 
presented. 

Chapter I of the thesis contains a brief introduction to 
the crystal spectra of rare earth ions and a summary of the 
proposed work on LiYF 4 :Eu^ + system. The over-all approach 
adopted to interpret the observed data has been briefly out- 
lined. Chapter II contains the theoretical background rel- 
evant to the present problem. Information on the crystal 
structure of LiYF^ and experimental detail relevant to the 
work are given in Chapter III. The experimental results are 
tabulated in this chapter. A detailed analysis of the expert 
mental data is presented in chapter IV. A summary of the 
results and their significance is discussed in the fifth and 
the final chapter. 

Some comments on the computer programming and calculation 
procedures are given in the appendices. 

I 



CHAPTER I 


INTRODUCTION 


Lanthanides form a group of chemically similar elements, 
which have a partially filled 4f shell in common. Neutral 
atoms of lanthanide group have Xenon like core of 54 elec- 
trons (Is 2 , 2s 2 2p 6 , 3s 2 3p 6 3d 10 , 4s 2 4p 6 4d 10 , 5s 2 5p 6 ), 

'N 1 4f electrons and two or three loosely bound outer 
electrons. The *3+' oxidation state for rare earths is 
the most stable, although the '2+' oxidation state is also 
encountered for some of them. In the optical spectra of 
rare earth doped crystals, one is primarily concerned with 
ionized rare earths, having only 4f^ electrons besides the 
Xe core. 

The existence of sharp bands in the absorption spectra 
of crystalline lanthanides was known from the work of 
Becquerel [l]. Further, they were known to show considerable 
Zeeman splittings [2]. However, it was Be the [3], who first 
proposed a highly successful crystal field theory to account 
for these observations. The absorption data of the aqueous 
solution of all rare earths (except Pm + ) were obtained as 
early as 1934 [4]- The interpretation of the complex spectra 



2 


of rare earths became possible only after Racah developed 
his powerful tensor operator techniques [5]. He introduced 
the concepts of tensor-operators, recoupling coefficients, 
seniority and coefficients of fractional parentage. Stevens 
[6] developed the equivalent operator technique for the 
theoretical interpretation of the spectra of complex systems. 
Satten for the first time applied these techniques to the 
4 f^ spectra of Nd^” in NdCBrO^) .9H2O. He, however, used 
Russell Saunders coupling [7]. It was subsequently realized 
that Intermediate coupling is more appropriate for these 
systems [ 8 , 9 ]. Since then considerable progress has been 
made in the interpretation of the spectra of lanthanide 
spectra [10,11,12]. 

The 4 f electrons, responsible for the spectroscopic 
properties of the lanthanides, are not the outermost elec- 

r\ /L 

trons. They are shielded by 5 s , 5 p electrons with larger 
spatial extension. As a result they interact rather weakly 
with the surrounding ions. This fact is responsible for 
the 'atomic like' spectra observed for rare-earth ions 
embedded in solids. The spectra of rare-earth ions in 
crystalline environment consist essentially of well separated 
groups of sharp lines. Each such group corresponds to the 
crystal field-split structure of the so called 'free ion 
levels'. Free ion levels of rare earths do not change much 
from one host to another, indicating the localized nature 



of the 4f electrons. The observed spectra of trivalent 
lanthanides originate from the transitions within the 4f^ 
manifold. 
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Spectroscopic study of rare earths is an interesting 
topic of research in itself. Lanthanides have proved to be 
good probes for the study of the crystalline field. Moreover, 
the use of rare earth doped materials [13] as gain media in 
lasers, has made this type of study more relevant. 

Tripositive europium (Eu ^ with six equivalent elec- 
trons in its 4f shell, is unique in the lanthanide series 

7 

in the sense that it has the ground term F q with total 
angular momentum J = 0. This imposes special restrictions 
on transitions involving the ground term. Further, the 
nuclear magnetic moment ofEu 3+ is quenched to such an extent 
that no NMR signal has been observed in any host by conven- 
tional NMR methods [ 14 ] . Laser action has been observed in 
Eu 3+ doped Y 2 0 3 and YV0 4 crystals [15,16,13]. 

Laser action has been observed in the LiYF 4 host when 

doped with Ce 3 *, Pr 3+ , Nd 3+ , Tb 3+ , Ho 3+ , Er 3+ , Tm 3+ etc 

34* 

[13], In LiYF 4 , rare earth ions substitute for Y with 
site symmetry D 2c j somewhat distorted to S 4 . There is strong 
evidence for the appearance of several magnetic dipole 
transitions in this system. This is somewhat unusual 
because magnetic dipole transitions, generally weak, are 
seldom observed. 
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In this thesis, we have undertaken the task of an in- 

depth study of the absorption and fluorescence spectra of 
3+ 

Eu :LiYF^. We have basically followed the procedure 

described by Wybourne [ll] and Hufner [13]. We obtained 

polarized absorption and fluorescence spectra at liquid 

nitrogen temperature in the visible spectral region. The 

absorption spectrum was obtained using a 1000W tungsten- 

halogen lamp, while for fluorescence studies, a Coherent 

Ar + laser and a Coherent ring dye laser with R 6G dye were 

used. Further, excitation spectrum coverincr the spectral 

range of R 6G dye was recorded to ascertain the position of 
5 3+ 

the D q level of Eu . The results have been analyzed in 
the framework of intermediate coupling scheme and crystal 
field theory. 

Weak interaction of rare earth ions with surroundings 
permits one to treat the crystal field effects as a small 
perturbation over the free ion Hamiltonian. In the theore- 
tical treatment of free ions, the non-central electrostatic, 
the spin-orbit and the configuration interactions are 
treated as perturbations to the central field Hamiltonian. 
Matrices, consisting of the sum of the matrix elements of 
the three interactions, calculated in LS basis states, are 
diagonalized to obtain the free ion energies and the inter- 
mediate coupling wave functions. The basis states |WUSLJ>, 
used to calculate individual matrix elements involve composite 
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orbital (L), spin (S) and total (J) angular momentum quan- 
tum numbers. U and W are additional quantum numbers intro- 
duced by Racah for a complete classification of the states 
N 

of the 4f configurations. The angular parts of the matrix 
elements can be calculated exactly, while radial parts are 
treated as adjustable parameters in an effort to match 
experimental and calculated energy positions of ’free ion* 
levels . 

The free ion levels are (2J+1) fold deqeneratei This 

degeneracy is wholly or partially lifted by the crystal 

field depending upon the symmetry of the field. Matrix 

elements of the crystal field are calculated using free ion 

wave functions as zero order wave functions. Basis states 

chosen are the linear combinations ja JJz>, denoted by 

|a JJz n> where a stands for the rest of the quantum 

'th 

i is the i x irreducible representation of 
the point group under consideration. Again, radial parts 
of the matrix elements are treated as adjustable parameters. 
We adopted this approach and used computer programmes to 
minimize the r.m.s. deviations between the theoretical and 
experimental positions of the energy levels and obtained the 
free ion and crystal field parameters for the system 
Eu^rLiYF^ in the intermediate coupling scheme. Subsequently 
the J-mixing terms were included in 'the calculations. 
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CHAPTER II 


THEORY OF CRYSTAL SPECTRA OF RARE EARTH IONS 


1. INTRODUCTION 

Trivalent lanthanides possess electronic configura- 
tions of the type (Pd core) 4 ^ 1 4f^, 5s 2 5p^, where 1<N<13. 
With the increase in the atomic number within the 
lanthanide series, the nuclear charge and the number of 
electrons in 4f shell increase. Due to imperfect shielding 
[l,2] among 4f electrons, this leads to an increased nuclear 
attraction for the f-electrons. This potential well near 
the nucleus reduces the spatial extension of the 4f wave- 
function to such an extent that its peak lies well 
within the spatial extension of the 5s and 5p wavefunctions . 
This reduction in the spatial extension of the 4f wave- 
function is named ' lanthanide contraction ' [2] . Figure (2.1) 
gives the radial distribution functions of the 4f, 5s, 5p 
and 6s electrons for Gd + as obtained from Hartree-Fock 
calculations of Freeman and Watson [3]. Thus 4f electrons 
are substantially screened from outside influences by the 
6s 2 and 5p^ electrons. This screening effect makes 4f 
electrons less sensitive to ligands in molecular and solid 
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FIG. 2.1 : RADIAL CHARGE DISTRIBUTION P^(r) AS A FUNCTION OF 
r FOR THE 4f, 5s, 5p and 6s ORBITALS [ll]. 
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complexes. In fact the spin-orbit and the residual 
Coulomb interactions for trivalent lanthanides are usually 
much larger than the ion-ligand interactions. These 
circumstances allow the crystal field to be treated as a 
perturbation over the ’free ion' interactions. 


2. FREE- ION HAMILTONIAN 


With the assumption that all the electronic shells 
except the 4f shell are complete and therefore do not 
contribute to the relative positions of the 4f energy 
levels, one can write the free ion Hamiltonian involving 
only the 4f electrons. 


H = -25 


T5 2 N 
n E 


v? - 


i= 1 


N 2 N 2 N 

E 2* e + E — + E ^( r .)(s. 

i=l r i i< j r i j i=l i i i 


Here Z* e is the screened nuclear charge, r^ is the 
coordinate of the 1 th electron w.r.t. the nucleus, r^ 

j. y. 4* h 

is the distance between i and j electrons and 
is the spin-orbit coupling function defined as 


S, ( r i> = • 

where U( is 
i th electron, 
meaning . 


~h 2 dU(rj) 

>m 2 c 2 dr i 

the spherical potential experienced by the 
Remaining quantities have their usual 


The first term in the Hamiltonian 1 represents the 
kinetic energy of the electrons, the second represents their 
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Coulomb interaction with the nucleus, and the third and 
fourth terms respectively represent the mutual Coulomb 
repulsion among the electrons and the spin-orbit inter- 
action. The last two terms are responsible for the energy 
level structure of the 4f^ configuration. 

It is impossible to solve the Schrodinger equation 
exactly for atoms containing more than one electron, 
analytically or numerically. Therefore one has to use 
some approximate models. The most commonly used approxi- 
mation for solving the Schrodinger equation for complex 
atoms is the central field approximation [ 4 ]. 

2.1 CENTRAL FIELD APPROXIMATION 

In the central field approximation, it is assumed 
that individual electrons of the atom move independently 
in a spherically symmetric central field (- created 

by the remaining electrons and the nucleus. Under the 
above approximation the Hamiltonian 1 takes the form, 

N -*2 2 N 2 

H = i=i c "^ v i + u(r i )] + hf - utr i> 

N 2 N 

+ Z ~] + Z S(r,) (s f .i,.) 
i<j r ij i=l n 1 1 ^ 

- H C F + H Coul + H ’^, 3 

The non-central residual Coulomb interaction H£ oul and the 
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spin-orbit interaction H^_q are smaller than the central 
field H c p by at least a few orders of magnitude and in the 
theory of atomic spectra they are treated as perturbations 
over the central field Hamiltonian. 

The time independent Schrodinger equation for the 
central field Hamiltonian 

H CF 'V CF = E CF V CF 4 

can be solved in the single electron approximation by 

writing the wavefunction ^CF anc * ener 9Y the form 



N 

E r _ = £ E. 5B 

Lh i=l 1 

Here is the spin wavefunction. k and a represent the set 

of quantum numbers (ni/n^) and (n JLm^m s ). (j)(k) for each 
electron satisfies 

[-^ V 2 + U(r)J |(k) = E <ji( k ) 6 

The equation 6 is similar to the Schrodinger equation for 
the hydrogen atom except for the replacement of the Coulomb 
potential -e 2 /r by the central field potential U(r). This 
can be solved by the introduction of polar coordinates 
(r, 9, (j)) and separating the one electron eigenfunctions 
into their radial and angular parts. The normalized bound 
state solutions for equation 6 can be written in the form 

<j)(k) = r \ R nt (r). Y t (0, $) 


7 
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where R n ^(r) is the radial function and depends on the 
central field potential energy function U(r). R (r) 

n *V 

satisfies the equation 


i- a- (r 
2 dr ^ 
r 


d R„ . (r) 


^_L) -Mmi %lM+ 2ffi [E-u(r)] 


R b (r) = 0 
nl 
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whereas spherical harmonics Y^ L (6, (f)) are defined by 


Y^(©> <t>) 


(-D mi [iakii fkii^ji/2 A (cos G) 

*e im ^ 
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and 



(x) 


9 m t/ 2 

(1-x 2 ) 

2 l .& I 


m<+l 


dx 




(x 2 -!) 1 
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The angular part of the one electron eigenfunction 
(Eqn.7) is identical to that of the hydrogenic eigen- 
function and can be evaluated exactly. Radial part can be 
evaluated only approximately in the absence of the knowledge 
of the potential function U(r^), Therefore it is customary 
to treat radial parts as adjustable parameters to fit the 
experimental data. 

The normalized antisymmetric solution of central 
field can be written in the de terminantal form as 
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fiU 1 ) 

^(cc 1 ) 

i N (^) 

^(oc 2 ) 

I 

3> 2 (a 2 ) 

— <M" 2 ) 

f 

$ 2 (« N ) 

- - - y« N ) 
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These wavefunctions will serve as zero order wave- 
functions. The eigen values of the central field wave 
equation are degenerate in quantum numbers m, and m and can be 
characterized by a sequence of quantum numbers (n^ 

(n 2 i 2 ) (n N l N ) which defines an electronic configura- 

tion. 

Once the zeroth-order wavefunctions are known one can 
use them as basis states to calculate matrix elements of 
residual electrostatic (Hq ou ^) and spin-orbit (H£_q) perturba- 
tions. The closed shells do not affect the relative splitt- 
ings in a configuration, but shift the whole configuration. 

So the optical spectra of rare-earth ions can be explained 
by considering electrons in 4f shell only. Even with this 
simplification, the direct method of calculating the effect of 
the perturbation (H ' = H£ 0ul + H s-0^ usin 9 determinantal 
product states as basis functions is tedious, except for 
very simple configurations. For example the configuration 
4f (Eu ), in which we are interested, has 295 free ion 
terms with total degeneracy 3003. This requires the evalua- 
tion of 295x295 = 87025 matrix elements and diagonalization 
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of a secular determinant of order 295x295. Moreover this 
method does not provide any physical insight into the 
structure of the configuration. The calculations are 
simplified if group theoretical and tensorial techniques 
are used. 

2.2 CLASSIFICATION OF THE STATES 

In the LS coupling scheme L and S define a Term. In 
complex atomic configurations, one frequently encounters 
more than one Term with same L and S values. Additional 
quantum numbers are then needed to label such Terms uniquely. 
Racah [5,6] developed a group theoretical approach to 
classify the states of ,£, N type configurations systematically, 
his method the complete set of states of the f^ configuration 
can be classified by their properties under certain groups 
of transformations. The irreducible representations of the 
group can then serve as suitable labels to classify the 
states. Complete classification of f^ configurations can 
be achieved using the unitary transformation group U 7 and 
its sub groups R 7 » R2 and such that 

u 7 dr 7 dr 3 dg 2 

The infinitesimal operators of the seven dimensional 

(k) 

unitary group U 7 are the tensor operators V q where 

-k < q < k and 0<k<6. V q k ^ commute with spin s. The 

N 

totality of the states of the configuration £ forms a 
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basis for a single irreducible representation of ^2^1* 

Since the representation does not depend on , the labeling 


of states by specifying (2J L+l) integers ^2 * * * * 

is equivalent to specifying the total spin S. The three 


dimensional rotation group R^ (a sub group of U^) has 

(k) 

tensor operators of odd rank as its elements. Its 
irreducible representation transforms as total orbital 
angular momentum L. This implies that the labeling of states 


with is equivalent to specifying L. Thus labeling by A 
and is equivalent to Russell-Saunders labeling. The L-S 
labeling is preferred over A anc * as it has more direct 
physical meaning.' 


Further classification is obtained according to the 

seven dimensional rotation group R^ and its sub group 

The irreducible representations of R^ are labelled by 

three integers (w^ to^) = W, such that 2 > to 1 9 anc * 

that of G 2 by two integers (u^ U 2 ) = U such that 

Using these labels as additional quantum numbers, the 

states of f N configuration can be uniquely defined with the 

exception of a few of the states with U = (31) and U = (40), 

where some doubly occurring LS states remain unseparated* 

These unseparated states are distinguished arbitrarily by 

giving an additional label X to these pairs of states. 

N 

Thus the states of configuration f can be completely 
specified by writing them as jf^XWUSL>. 
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Racah also introduced the concept of seniority numbers 

[6] which give an alternate way to classify the states of 
N 

f configurations and are convenient for the calculation of 
matrix elements of certain operators. In this approach the 
terms |SLJMj> of a configuration can be divided into two 
types. The first type of terms can be obtained by the 
addition of a pair of electrons in the | configuration, 
the other type can not be obtained in this manner. The 
second type of terms can be considered to appear for the 
first time. A unique number • v* called seniority number, 
can be assigned to the term appearing for the first time 
in a configuration, starting from the configuration 
The term S=0, L=0 will have v = 0. Classification of terms 
using seniority number is equivalent to specifying *W'. Thus 

i N 

the states denoted as |f X WUSL> can be equivalently written 
as | f ^ xvUSL>«. The complete classification of the terms of 
the 4f^ configuration is given in appendix-I. 

2.3 CHOICE OF BASIS STATES 

The matrix elements of the perturbation Hamiltonian 
(equation 3) can be calculated by choosing a complete set 
of basis states in any well defined coupling scheme. The 
Russell-Saunders (L-S) coupling is applicable when H£ QUi » 

H£_q and the j-j coupling scheme works when H c ou i <<H s-0* 

For the rare earth ions H s_ol H c 0 ul and therefore neither 
of the coupling schemes are appropriate. The coupling 



scheme commonly used for calculating the energy levels of 
rare earth ions is called the intermediate coupling scheme. 
The matrix elements of the perturbation Hamiltonian are 
usually calculated using the L-S basis states and the 
secular equation of the composite Hamiltonian +Hg__Q) 

is then solved. Since H^ ou ^ is diagonal in |SL> representa- 
tion and independent of J and J z » while H£_q is diagonal in 
J and independent of J z , we calculate the matrix elements 
with | SLJ> bases. This choice breaks- up the secular equa- 
tion into smaller secular equations, one for each J value. 

Now the secular equation for a specified J takes the form 


<1 |H r j 1>-X 

<2 1 H ' 1 1> 

< i | H ' 1 1> 

f 

I 

<n|H ' | 1> 


<1 1 H ’ l 2> . 

< 2 | H 1 I 2>"X» 

< i | H * | 2> . 


<1 |H ' J n> 
<2 | H 1 1 n> 
< i | H 1 |n> 


<n | H ' 1 2> 


.<n|H' 


where <i| = <S^L^j| and |j> = 
and H’ = H£ oul + 

The roots of this equation € 2 6 n 9 ive the 

required energy correction. Once the secular equations are 
solved eigen functions j can be determined. 


The problem now reduces to the calculation of the 
matrix elements <S i J J z j H£ oul + H£_q| J J 2 > . 

Since both the perturbation Hamiltonians H£ qu1 and H£_q 
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are the functions of the coordinates of individual electrons 
and the states jS L J J z > are of composite nature, this require 
the calculation of the matrix elements to be carried out 
in de terminantal product states. This will be a difficult 
task for complex configurations. Tensor operator techniques 
developed by Racah [6,7] make these calculations practically 
simpler. Using these techniques calculations are done 
directly in composite states by an appropriate transforma- 
tion of the Hamiltonian and using the concept of fractional 
parentage coefficients [6]. 

2.4 TENSOR OPERATORS [6] 

(k ) 

Spherical tensor operator T' of rank k is a set of 
(2k+l) operators Tq k ^ (q = -k, -k+l,. . .k) having trans- 
formational properties similar to those of spherical 
fk} (k) 

harmonics Y^ ' • This means T* 1 will satisfy the same 

(k ) 

commutation relations as those of Y' with the total 

M 

angular momentum J. 

[(J x + J y )» T q k) 1 = L(k+ q)(k+q+l)] l/2 .....13 

and 

[J z , T^ k) ] = q T^ k) 

In the JJ 2 representation, the dependence of the matrix 
elements of tensor operators can be obtained by Wigner- 
Eckart theorem. 

<Y JJ z l T q k ^ I = (-D J - J Z <-TJ I |T (k) |.|yj»> 


J 


( 


-J 


z 


k 

q 


j' 

j* 

z 


) 


14 
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The quantity, <y J||T (k) ||y J'>, which ls independent of 

^z* an< ^ q » called the reduced matrix element, 

( J k J * \ . 

x „ T , ' 1S a 3 “ J symbol and y stands for the rest of 

z q J z 

the quantum numbers on which T^' k ^ does not act. 

The reduced matrix elements <y J| | | | y' J‘> can be 

evaluated if the J values and the nature' of are known. 


Another tensor operator, called 'product tensor 

operator' is encountered in the calculation of the matrix 

elements of H£ and H| . A spherical product tensor 
(k) . 

operator Xq is defined as the product of two spherical 

tensor operators and 

q l ^2 

4 KL * T< k l> U< k 2> (2K+1) (-l) k l“ k 2 +Q ( k l 

q l q 2 1 2 q x 

- ^T (k l). U fk2) ] ™ 



satisfies the same commutation relations as ' S 

Q . 

If kj=k 2 =k, the scalar product of two tensors can be 
defined as 



(?(k)^ u(k)j = 


I (-l) q T^ k ^ l/ k ^ 
q q -q 


or 

(T (k) .U (k) ) = (-l) k (2k+l) l/2 ^T (k ).U (k) j° . 

Again Q dependence of the matrix elements of Xg k ^ 
be given by Wigner-Eckart theorem. 


16 


can 
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(Xj 1 j 2 JJ 2 |x^ir j- j- J’J’> 


= (-l) J “ J z ( J )(tj 1 j 2 J||X (K) ilx’ ji JAJ’) 17 

-J Q J' ^ 

Z Z 

The reduced matrix element can be expressed in terms of 

the product of two reduced matrix elements in the followinq 

(ki) (k 2 ) 

way (If T and U operate on two different parts 1 and 


2 of the system respectively). 
<Zhh J MX (K) | Ix’JiJA J> = T 


(ki) 


J> = I I'd' ji> 


(k 0 ) 


<x" 1 2 I ^ U I ix'i’ > -"f[(2J+l)(2K+l)(2J' + l)]>j 2 y 2 k 


h Ji k : 


J J’ K 


where the last factor is a 9-J symbol. 


For the case k-^ = ^2 = k, K = 0. Using the definition of 
scalar product of two spherical tensors we have 

<x jp 2 ji jjj'jp 


JV 

z z 


X L„ < T | T (k) | lx" Ji>< z" j 2 l lu (k) I lx; jp 

w 

where the factor in curly brackets is a 6-J symbol, 
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The 3-J and 6-J symbols are tabulated by Rotenburg et al [8]. 


2.5 MATRIX ELEMENTS FOR COULOMB INTERACTION 

In equation 3 for the free ion Hamiltonian, we have defined 
the non-central part of the Coulomb interaction as 
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HA..., = e [- 


Coul 


i=l 


, N e 2 
- U(r ± ) ] + E 


.20 


i<j i 3 


The terms in the bracket are purely radial and centrally 
symmetric. They shift the configuration as a whole and do 
not participate in the relative splittings of the terms in 
a given configuration. They can be ignored for the present 
problem. The term E e 2 /r i j removes the degeneracy of the 

i<3 J 

terms partially. The tensor operators corresponding to the 
last term can be obtained by transforming it in such a way 
as to separate the radial and angular variables. The factor 
l/r^j can be expanded in a series of Legendre polynomials 


ij 


= — r r — = (rf + r? - 2 r_- .r-j Cos to) 

I “ “I i J i J 

|r i " r j ! 

r, k 

= A -Tc TT P * < Cos 


- 1/2 
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Here r^ and r^ denote the smaller and greater of the 
magnitudes of the vectors r^ and r j , and to is the angle between 
the vectors r^ and r j . Using the theorem of addition of 
spherical functions, P k (w) can be written in terms of 
spherical harmonics Yq k ) (©i^) and Y^ 1 ^ (©j(J)j) . Let us also 
define a single particle tensor operator ^ (i) as, 


c (k) (i) _ ( J2-. 

u q ~ v 2k+l 


■) 1//2 Yi k ^ (i) 
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Thus the expression for e 2 /r^j becomes 


« 2 ^i 3 = e2 | 5KT TTZ1 * Y< k) (9i+i). 


.(k) 


q=-k 
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e 2 l 2 C< k) (i) C< k) (j) 

k r> +1 q=-k q q 


e 2 I 


r > 


**<A 


I (-l) q C* k) (i) C* k, (j) 


e 2 E 


r > 


( C (k) r (k) x 

k+1 ^ i * j ; 
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The matrix elements for the residual Coulomb interac- 
tion can be evaluated by using equation 19 for the scalar 
product of two tensor operators. 


Using | X SL> basis states (where *£ stands for the rest 
of the quantum numbers) the matrix elements can be written 
as 


<nf N TSL| i ^ J §^j|nf N x* S' L*> 


= <n f\s L| e 2 Z £ 


i<J k r> k+1 


(c{ k) 


= E f u <nf N X SL, nf N x* S' L'>. F k 
k K 


. Cj k ^ ) | n f N S»L*> 
(nf,nf) .....24 


Here angular parts f. of the matrix elements are defined as 


,25 


f k = <nf N .xSL| i E_ } C^. cj k Mnf N x' S' L*> 

= 6 SS , ^ ll ,(-D L ^ 3 ^(<3| |C (k) 1 |3>) 2 .. 

and Slater's radial integrals F are defined as 

F* = F k (nf ,nf ) = e 2 f £ ~^[[R nf (rj.) .R nf (r^ ) ] 2 dri drj ..26 
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The reduced matrix elements <3.| |c^ k ^ | ]3> can be calculated 
using the expression 

<l| |C (k) | |t'> = (-1) 1 XV(2t+l)(2i*+IT (* k J) 

Thus 

<3 1 |C^ | 1 3> = - Y8T” X (q o q) 27 

The triangular condition and the condition k+t+JL’ = 2g 
(g -integer) for non zero values of 3-J symbols restrict 
the values of k for non-zero matrix elements to 0,2,4 and 6. 


To avoid the occurrence of large fractional coef fi- 
le 

cients for F in the evaluation of matrix elements, Condon 

k k 

and Shortley redefined f and F as 


f,. = 


f k r, _ F k 


k " T3£ and F k = 


.28 


Where D. are the denominators listed in the Table 1 and 
k 

2^ of reference [4]. 

The operator (C^ k ^. C^ k ^) does not have the transfor- 
mational properties of groups Ry and G2, used to classify 
the states of the configuration, though it is a scalar with 
respect to the group R3. This makes the calculation of 
matrix elements difficult for systems having more than two 

equivalent f electrons. To overcome this problem Racah 

lc 

defined a new set of parameters e^ and E , which are linear 

tr 

combinations of the parameters f and F fc . The matrix elements 

in terms of Racah's parameters can be expressed as 

3 k 

E = E e. E K 

k=0 k 
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Where e^, the angular parts are related with the old 
parameters f^ by the expressions: 

e Q = f° = N(N-l) 

e x = (9/7) f° + (1/42) f 2 + (1/77) f 4 + (l/462) f 6 

e 2 = (143/42) f 2 - (130/77) f 4 + (35/462) f 6 

e 3 = (11/42) f 2 + (4/77) f 4 - (7/462) f 6 30 

and the new radial parameters E^ are related to FjJ s as 

• E° = F q - 10 F 2 - 33F 4 - 286 F 6 
E 1 = |(70F 2 + 231F 4 + 2002F 6 ) 

E 2 = ^(F 2 - 3F 4 + 7F 6 ) 

E 3 = 3(5F 2 + 6F 4 - 91F 6 ) 31 

Nelson and Koster [8], have calculated and tabulated the 

N 

angular parts e^ for all t configurations. The radial 
quantities E are treated as adjustable parameters in the 
interpretation of experimental results. 

2.6 SPIN-ORBIT INTERACTION 

The spin-orbit Hamiltonian defined in equation 3 as 
, N _ 

H S -0 = i S i ^(r i )(s i .« 4 ) 32 

is a simple example of the product tensor operator of 
two tensors of rank one; one operating only on the spin 

and the other operating only on the spatial coordinates. 

2 2 2 
HI ~ commutes with J and J but not with L and S . 

1 S— 0 z 


As a 
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result the matrix elements are independent of J and diagonal 
in J but mix the states of different L and S values. 


The matrix elements of Hg_Q in the basis functions 
|f^*CSLJJ z > can be obtained as 

<f N TSLJJ z |H£_ 0 jf N -£» S ' L ' J ' J z > 


2 Z 


X <f N TSL I £ (s i .%)if IN x'S'L , > 
i=l 


N 


J*. 


,33 


Here, the spin-orbit radial factor ^ is defined as 

^nf = ^°° R ; nf (r) *^ (r) dr 34 

o 

The last factor can be further evaluated as 
<f N TSL| E (si.ti) |f N x* S ' L,> = [3(3+1) (2x3+1) 3 1 / 2 

• <f N XSL| | V^ 11 |f N -c r S'L’> 

= V84.<f N X SL| IvO 1 )] |f N r»S , L , >..35 

The reduced matrix element <| |V^^j j >can be evaluated 
using coefficients of fractional parentage[6] ; the results 
can be written as 

<f N TSL| |V^| |f N x’S’L^ = 3N[(2Sfl)(2S»+l)(2L+l)(2L , +l)] 1 / 2 
x( _l)l/2 +3+S+L # E( vpa(f ) (VVy'f)(-l) S+L 

y 

r s i so 

* (1/2 S 1/2J * 

(y^l^) are coefficients of fractional parentage. The reduced 
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matrix elements of for all p^, and f^ confiaurations 

are calculated and tabulated by Nielson and Koster [fl]. The 
triangular conditions for the 6-J symbols restrict the non- 
zero matrix elements to those with 

As = 0 , +1 and AL = 0, +1 37 

2.7 CONFIGURATION INTER AC TION 

Free ion calculations carried out 'within a single confi- 
guration approximation usually leads to a disagreement (of 
the order of hundreds of wave numbers) between calculated and 
experimentally estimated free ion energy levels, even if the 
radial integrals are treated as parameters. To consider the 
complete configuration interaction, one has to take into 
account very large number of configurations including those 
in the continuum. This will require handling unmanageably 

large size matrices. Therefore the commonly adopted method 

* 

to treat this type of interaction is to incorporate the most 
dominant effects of the interaction within the states of the 
ground configuration. The most dominant type of configuration 
interaction is via the Coulomb interaction. For the JL N 
configurations, this effect has been the oreti cal ly studied 
by Rajnak and Wybourne [9] and also by Racah 
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and Stein [lO]. Following are the most likely configura- 
tions which can perturb an configuration. 

(1) a n_1 r 

(2) fc N - 2 £ 2 and f-2 £ . 

( 3 ) e ' 4 t ' +1 £ N+1 

(4) i' 4i ’ +1 £ N + 1 

(5) £ ' 41 ' {_ N+2 and £ ' 4t ’ +1 i" 4 *’’ +1 A N+2 

The last three configurations correspond to core 
excitations and their effect is to shift the centre of 
gravity of the perturbed configuration as a whole. They 
do not contribute to its s true ture, hence these can be 
neglected from the present discussion. 

Rajnak and Wybourne [9] showed that the configura- 
tion interaction through the Coulomb field adds a correc- 
tion term 

• • • • • 38 

^T7 

to each of the matrix elements of the perturbed configuration. 
Here m_* reoresents perturbing conf iouration states which 

are separated by an energy E^- from the perturbed configura- 
tion and G is an operator representing the Coulomb energy 

r\ 

of repulsion, I e /r. . . They were able to show by using 
i< j J 

tensor algebra and tedious summations, that the above term 
can be expressed in terms of the matrix elements of the 


<t I 

C(W) = T 
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tensor operators l/ k ^ 

C(H > 4 , ‘) = I P k (f N ^| E (U.^ .U? k ^ ) |f N U>» ) 39 

k=0 i<j 1 j 

k 

Here P contains radial integrals and excitation energies. 

For even k the angular part of the equation (39) is 
having exactly the same form as the coefficients of F k in 
the matrix elements of Hence treating the radial 

integrals of the matrix elements of as adjustable 

parameters automatically includes this part of the configu- 
ration interaction. 

The terms with k odd, are not taken into account even 
when the Slater integrals are treated as adjustable para- 
meters. These can be included only by introducing addi- 
tional parameters. Rajnak and Wybourne used properties of 
Casimirs operators [6] to replace these terms with an equivalent 
expression . 

6(HT) [a .L(L+l) + p G(G 2 ) + y G(R ? )] 40 

Where GiG^) and G(Rj) are the eigenvalues of Casimir's 
operators for the groups and R^ used to classify the 
states of the f^ configuration. These are tabulated in 
Tables 2.6 and 2.7 of reference (ll). 

a, P and y are associated parameters corresponding to 
a particular linear combination of the P k (odd) parameters 
and are to be treated as adjustable parameters in order to 
improve upon the agreement between theoretical and experimental 
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results . Wybourne [ 12] has generalized the configuration 
interaction to higher orders of perturbation, but this needs 
introduction of additional parameters and in the absence of 
sufficient experimental data, it is not possible to consider 
these types of effects. Rajnak and Wybourne [13] have also 
examined configuration interaction via spin-orbit interac- 
tion. This simply screens the spin-orbit coupling radial 
integral and is taken care of by parameterization of spin- 
orbit interaction. 

2.8 OTHER INTERACTIONS 

Though the normal Coulomb, spin-orbit and configuration 
interactions take care of the major contribution to the 
energy of the states of the configurations encountered in 
the rare earths, some other interactions are also present; 
for example, orbit-orbit, spin-spin and spin- other orbit 
etc. Most of these magnetic interactions do not make any 
substantial change in the energy contribution while their 
consideration increases the complexity of the calculations. 
It may however be pointed out that the parameterization 
technique already discussed does take care of some of these 


effects 
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3. CRYSTAL FIELD INTERACTION 

In an ionic crystal the electrons occupy orbitals that 
are highly localized about the ions in such a way that any 
electron can still be associated with a particular ion. 

Each electron, however, feels the influence of the electrons 
belonging to the other ions (a repulsion) and of the nuclei 
belonging to the other ions (an attraction). We can take 
this influence into account by considering that the electrons 
of the individual ions are subjected to the action of a 
'crystalline field'. 

The 4f shell of rare earths is incomplete and therefore 
has an aspherical charge distribution. The crystal field is 
completely external to the ion and has a definite symmetry, 
invariant under time reversal. The crystal field experienced 
by the ions distorts the closed shells of the RE ions and 
removes, to a certain extent, the (2J+l) degeneracy of the 
free ion levels. In the case of an even electron system the 
degeneracy can be completely lifted by an electric field 
(crystal field) of low enough symmetry. For an odd electron 
system, at least a two fold degeneracy remains, irrespective 
of the symmetry of the field. This is known as Kramers' 
degeneracy [ 13 ] . This degeneracy may be lifted only by 
magnetic interactions. It is possible to predict the extent 
to which the degeneracy of the energy levels can be lifted 
and also the transformation properties of the eigen functions, 



31 


if the symmetry of the crystal field is known. The actual 
evaluation of the crystal field splittings, however, requires 
the solution of the crystal field Hamiltonian. The 
Hamiltonian for an ion in a crystal can be written as 


H = H P + H ' 

F cry 
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Where Hp is the free ion Hamiltonian already considered and 
H* represents the interac- 

cry 

environment and is given as 


H^ry represents the interaction of the free ion with its 


e -* P( R) 


H* = - Z f -4 
cr ^ i I R - r. 


dX 


.... 42 


Here P( R) is charge distribution which interacts with a 

J.U 

given electron, is the position vector of the i 
electron with respect to the nucleus. 


The Hamiltonian H' , can be expanded in terms of the 

cry 

spherical harmonics in a fashion similar to the H£ ou ^ term 
of the free ion Hamiltonian 

Thus H iry " - l /-! dt 

i R -hi 


= P k (CosfR,?^) -%-j- dx 


> r. k 


= /f<R) 2fe- ( - 1)q Y q k)(e i- ( t | i ) ^rr dT 

r k 

= -l ej P(R)(-l) q c£ k) (®, 4 >) d T 

ikq 1 q l l q r > A 
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= X 8* 43 

kqi M M 

v 

Where the crystal field parameter B to be treated as an 

M 

adjustable parameter is defined as 

k 

= -e/(-l) q f(R) (©,({>) dt 44 

r> 

3.1 NUMBER OF TERMS OCCURRING IN THE EXPANSION OF H * 

1 ~ - - - . . - - ciy 

The allowed values of k and q in equation 43 depend on 
the site symmetry of the ion under consideration. The 
crystal potential must possess the symmetry of the lattice 
under consideration, that is, it must be invariant under the 
operation of the point group. In the case of LiYF^JEu • 
Eu^has site symmetry of S 4 « This means that the potential 
function in this case must be invariant under the operations of 
thegroup S 4 . The symmetry operation S 4 can be expressed by 
the transformation 
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Using the cartesian representation of spherical harmonics and 
the above transformation, the H^ ry . for 4f electrons can be 

obtained . 


H 


+ Z „ 
k=4,6 


(e ) - B° C (o) + B 2 C (2) + B 4 C (4) + B*> C<*> 
cry' S 4 " B o U o + B o u o o o o o 

[B 4 (C 4 ^ } + C<*>) + i B^cJ k) - C^)1 46 
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Here the restriction k<6 applies because we are considering 

4f electrons. The odd terms are neglected because their 

matrix elements within the states of the same parity are 

identically zero. In writing Eqn.46 we have made sure that 

H' is Hermitian. This requires 
cry 

= (-l) q B^* 47 

From the properties of spherical harmonics, operators 
(cj^ + are real and hence coefficients B^ will also 

be real. The operators (c| k ^ - C^ ) are imaginary, there- 
fore, to make H^ r y hermitian and B^ k real, * i * has been 
introduced explicitly. 

The first term in the expansion with k = q = 0 is 
spherically symmetric. It does not contribute to the 
splitting of the free ion terms and may, therefore, be omitted 
in the calculations. However, if the precise configuration 
interaction is to be considered, it can not be neglected. 

Once the H' has been obtained in tensorial form, its 
cry 

matrix elements can be calculated by selecting a suitable 
set of basis states. 

3.2 CHOICE OF BASIS STATES 

The point group involves a 90° rotation followed 
by a reflection in a plane perpendicular to the rotation 
axis. It is characterized by four irreducible representa- 
tion for the single group, two of which are degenerate. 
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The relevant character table of the point group S 4 is shown 
in Table 2. 1« In this table a(b) under the column time 
reversal indicates that the basis functions for the represen- 
tation in question are real(complex) . Tables 2*2 and 2. 3 list 
the relevant portion of the multiplication table and the 
full rotation compatibility table respectively. We note 
that the representations and are degenerate and for 
the purpose of calculating the crystal field splittings, we 
need to consider only one of them. 

H cry does not commu te with J and J , so it mixes the 

states with different J and J values and hence J and J 

Z Z 

are no more good quantum numbers. However, the matrix 

elements of the crystal field interaction, between wave- 

functions belonging to different irreducible representations 

are zero. Thus the irreducible representations of the point 

group act as good quantum numbers. It is always possible to 

i N 

find some linear combinations of | f J J^> states having 
transformation properties like those of the irreducible 
representations of the point group under consideration. 

These irreducible representations can be used to label the 
basis states |f^ J J I -7 . > where P. is the irreducible 
representation. The appropriate J 2 ~linear combinations can 
be obtained through the use of the projection operators 
Pj^l.14]. The projection operator for a group G, 

when operated on a set of functions |JJ Z >, 


selects those 
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TABLE 2.1 

CHARACTER TABLE FOR GROUP [ 15] 


S 4 

E 

E 

S” 1 

b 4 

r=x 

b 4 

C 2 

=2 

S 4 

S 4 

Time 

Inver- 

sion 

Bases 

E 

1 

1 

1 

1 

1 

1 

1 

1 

a 

S z 

r 2 

1 

1 

-1 

-1 

1 

1 

-1 

-1 

a 

Z or XY 

r 3 

1 

1 

i 

i 

-1 

-1 

-i 

-i 

b 

-(s x + iS y ) 

r 4 

1 

1 

-i 

-i 

-1 

-1 

i 

i 

b 

or 

i (X - iY) 

(S x - iSy) 


or 


-i(X + iY) 


TABLE 2.-2 

MULTIPLICATION TABLE FOR S 1 [ 15] 
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TABLE 2,3 

FULL ROTATION GROUP COMPATIBILITY TABLE FOR THE GROUP 
S 1 s CONFIGURATION 4f 6 fl5l 


n 


D I 

n 


03 + n > 


n 

+ 2 P 2 

+ <r 3 + 5 ) 

D 3 

n 

+ 2 fJ 

+ 2 (T 3 +IJ) 

°4 

3 ^ 

+ 

+ 2 (P 3 + (J) 

D 5 

3 Pi 

+ 2 P 2 

+ 3(T 3 +I 7 ) 

°6 

3 Pi 

+ 4P 2 

+ 3(G 3 + Q) 

D 7 

3 Pi 

+ ^2 

+4<r 3 + [p 

D 8 

5tl 

+ 

+ 4 (P3 + IP 

D+ 9 

5 Pi 

+ 4\7 

+ 5(P 3 + 1 J) 

D to 

S E 

+ 6 

+ 5(T 3 +G) 

D L 

5 Pi 

+ 6q 

+5<r 3 +tj) 

CM 

+ cT 

7 Pi 

+ 6 ^ 

+5(r 3 + 14 ) 



37 


*fc h, 

functions which transform according to the k row of the 
i^" 1 irreducible representation of the group G. For a point 
group of type S 2n 

|cr J J z > = 1 cr J J z >. 6(J Z » p(i»k) +nu + 2np) 48 

Where a = 1 for even parity configurations 
= 0 for odd parity configurations 
p = 0 , +1» +2» +3* • «•• 

For S 4 , n .= 2 and 

k = 1 for all irreducible representations of S 4 « 
p(i,l) = 0 and 2 for i = 1 and 2 respectively 
= + 1 for i = 3 and 4. 

Table 2.4 shows the basis functions for the irreducible 
representations of S 4 for the 4f 6 configuration. Numbers 
under the column N(fi) represent the number of times the i 

+ if * § 

representation appears when Dj is broken into 
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TABLE 2.4 

THE J^-BASES FOR THE IRREDUCIBLE REPRESENTATIONS OF THE POINT 
GROUP S 1 APPLICABLE TO Eu 3 * 


J 

n 


1 




No. of 
Stark 
Compo- 
nents 

£ N(R) 
i 1 

J z 

Nirp 

A 

nTT|7 

A 

imp 

0 

0 

1 





1 

1 

0 

1 



1 

1 

2 

2 

0 

1 

+2 

2 

1 

1 

4 

3 

0 

1 

±2 

2 

1,-3 

2 

5 

4 

0 ,+4 

3 

±2 

2 

CO 

I 

•o 

2 

7 

5 

0,±4 

3 

± 2 

2 

0,-3, 5 

3 

8 

6 

0,+4 

3 

+2 ,+6 

4 

1,-3, 5 

3 

10 

7 

0, + 4 

3 

+2, +6 

4 

1,-3, 5,-7 

4 

11 

8 

0 ,+4,+8 

5 

+2 ,+6 

4 

1,-3, 5, -7 

4 

13 

9 

0 , +4 , + 8 

5 

+2, +6 

4 

1,-3, 5, -7, 9 

5 

14 

10 

0,+4,±8 

5 

+2, +6, +10 

6 

1,-3, 5, -7, 9 

5 

16 

11 

0,+4, + 8 

5 

+2, +6 ,+10 

6 

1,-3, 5, -7, 9, -11 

6 

17 

12 

0 ,+4 ,+12 

7 

+2, +6 ,+10 

6 

1,-3, 5,-7, 9,-11 

6 

19 
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3.3 MATRIX ELEMENTS FOR H' 

■■ cry 

Matrix elements of can be calculated by tensorial 

methods. A typical matrix element will be given by 

<f N SLJJ |H* lf N S*L'J»J’> = E B?f<f N SLJJ |ul k ^|f N SL’J’J^ 

ZLXy Z ICCf^ Z Cf Z 

x <f I |c (k) | 1 f > ....49 

Further J z dependence of the matrix element of the tensor 
operator can be obtained using the Wigne r-Eckart theorem. 

N | (k)i N / ^ 

<f w SLJJ z |U^ k; |f 1N SL ' J ' J^> = (-1) z V ) x 

~ J z q J z 

<f N SLJ | | k ^ | | f N SL'J*> 50 

Where 

<f N S L J 1 | U ^ k ^ 1 1 f N SL' J’> = (-1) S+L ' +J+K V‘[(2J+l)(2J'+l)] 

j . , 

LL’ L S 

The doubly reduced matrix elements <f N SL|jl/ k ^||f N SL'> 
are calculated and tabulated by Nielson and Koster [8]. 
These expressions can be used to calculate matrix elements 
of the crystal field Hamiltonian within a multiple t or 
between multiplets. 

In equation 49, the terms with q = 0, (tT k ^) are axial 
“t 0 rms and cfiva th© non zaro matrix alamants for a Qivan 




<f N SL i lu< k >| |f N SL'>... .51 
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whereas the terms mix states for which J' - J = q. 

Because of this mixing, J and are no more good guantum 
numbers. The odd terms in do not contribute to the 

crystal field splittings. They are, however, responsible 
for the appearance of electric dipole transitions in the 
spectra of trivalent rare earth ions by mixing configurations 
of opposite parity. 

4. SELECTION RULES 

4.1 FREE ION SE LECTION RULES 

The electric dipole transitions among the levels of a 
N 

4f configuration are forbidden by the parity rule. The 
magnetic dipole transitions though allowed, are usually 
weak. For an ion placed in a crystal lacking a centre of 
inversion, mixing of configurations of opposite parity can 
take place due to the presence of odd terms of the crystal 
field Hamiltonian. The resulting electric dipole transitions 
are called 'forced electric dipole transitions'. Judd [ 16 ] 
and Ofelt [l7] have indeoendently developed the theory of 
forced electric dipole transitions. They have obtained the 
following selection rules for the magnetic dipole and for the 
forced electric dipole transitions. 

For electric dipole transitions: A J<6 

For transitions involving J=0 state, the following additional 
selection rules also hold. 
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(i) J = 0 J’ = 0 forbidden 

(ii) J = 0 J’ odd weak 

(iii) J = 0 J’ = 2,4,6 allowed 

For magnetic dipole transitions 

(i) A J = 0, +1 

(ii) J = 0 J* = 0 forbidden 

J - mixing by the crystal field can relax these selection 
rules . 

4.2 CRYSTAL FIELD SELECTION RULES 

The selection rules mentioned in the earlier section 
need to be supplemented due to the additional constraints 

imposed by site symmetry considerations. The character table 
for the ooint group S 4 (Table 2.l) shows that the oper- 
ator z and (x+iy) transform according to fg and re- 

presentations respectively . An electric dipole transition 
between levels characterized by the irreducible representa- 
tions ^ and fj will be allowed in it(cx) polarization 
provided the product fi p 2 ( [3,4) fj contains the symmetric 
representation fl- For example the product h r 2 (i contain 

Pi and therefore (2 transitions is allowed in it- 

polarization but the product V 1 V 2 V 2 does not contain and 
hence fi— - (2 is not allowed in cr-polarization. Using these 
considerations, the E.D. and M.D. rules listed in Table 2.5 


can be obtained. 
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TABLE 2.5 

FORCED ELECTRIC DIPOLE AND MAGNETIC DIPOLE SELECTION RUL ES 

FOR LiYF.:Eu 3+ 

4 — 


E.D. 

n 

R 

R.4 

a 

— 

71 

(5 

r 2 

1 

— 

a 

"7 

3,4 

a 

a 

71 


M.D. 

Fl 


[ 3,4 

tl 

a 

- 

71 

fl 

- 

a 

71 

f3,4 

71 

7E 

O 


The Pg and representations 


are degenerate 
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CHAPTER III 


EXPERIMENTAL DETAILS 


1. GENERAL 

This thesis deals with the observation and inter- 
pretation of the absorption and fluorescence spectra of 
Eu 34 doped in LiYF 4 single crystal. The single crystal 
used in this work was obtained from Dr L.E. Erickson of 
the National Research Council, Ottawa (Canada). The 
dimensions of the oriented and polished sample were 3.8 mmx 
3.1 mm x 2.1 mm with the c axis parallel to the 3.1 mm edge. 
Eu 34 concentration in the sample is approximately l‘/» 

2. CRYSTAL STRUCTURE [l] 

LiYF^ belongs to the tetragonal scheelite (CaWO^) type 
family of crystals with 

a = b 4 c 
a = p = y = 90° 

6 > 

The scheeli tes have space group C 4h (i^/a). In the 
isomorphism between LiYF 4 and CaW0 4 , Li 4 corresponds to W 6+ , 
Y^ 4 corresponds to Ca' 34 and F to 0 . The unit cell 

contains four molecules. The positions of the atoms in the 
unit cell are [2] 
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L : (0,0,0); (0,l/ 2 ,l/ 4 ) ; Cl/2.1/ 2>1/2) ; 

Y : (0,0,l/2);(l/2, 0,1/4); (l/2,l/ 2f0) . ( 0tl / 2f3 / 4 ) 
F : (x,y,z); (x,y,z); (y,x,z); (y,x f -j 


(x,y+ 1/2,1/ 4-z) ; (x, l/2-y, 

(y , 1 / 2 -x , z+ 1 / 4 ) + bcc. 


1 ^ 4 '“ 2 ) ; (y, x+l/2, z+l/4) ; 


Where x = 0.25, y = 0.15 and z = 0.075 n 

Crystallographic 

information on some of the scheelites i, 

-‘-is ted in Table 3.1. 

Figure 3.1 shows the positions of Li and v • 

- ions in the unit 

cell of LiYF A . The fluorine ions have n n + . 

^ x been shown in 

the figure. If we consider the site (i/o , , 

vj - / 2 »1/2,1/2) in 

Figure 3.1, we notice that the Yttrium i , , 

0n there is surrounded 

by 8 Li ions: 4 at the corner of a square 
1 re with a Y-Li 

distance of 3.72 A°, and 4 others at + 

Lne corners of a 

tetrahedron with a Y-Li distance of 3.70 ,0 

A . These are 

marked 1 to 8 in Figure 3.1. The f 0 „_ 

OUr fluorine ions 

surrounding a Li site lie on the vertices ^ . . 

s °f a somewhat 

distorted tetrahedron. One fluorine i 0 n -f~ „ j., 

xrom each of the 

S'CLiF.) tetrahedra contribute to the «,* , 

4 e ight fluorine 

neighbours of a Y ion, which is nearly a dnn a . 0 0 

uae cahedron. Fig 3.z 

shows the projected positions of the fi.,-.. 

u orme neighbours of 

a Y ion. The position coordinates of tho . . 

ne first nearest neigh- 
bours and the second nearest neighbours 

are listed in 

Table 3.2. Table 3.3 gives the radii of Li + y 3+ Eu 3+ and 
F*~ ions . 
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TABLE 3.1 

LATTICE PARAMET ERS OF LiYF 4 AND SOME OTHER SCHEELITES [3] 


S.No. 

Crystal 

C(A°) 

a (A 0 ) 


c/a 

1. 

LiYF 4 * 

10.94 

5.26 


2.08 

2. 

LiYF 4 ** 

10.85 

5.16 


2.10 

3o 

CaMoO„ 

4 

11.194 

5.1554 


2.171 

4. 

CaW0 4 

11.376 

5.243 


2.170 

5. 

PbMo0 4 

12.0165 

5.0165 


2.121 



TABLE 3. 

2 



POSITION COORDINATES OF FLUORINE IONS IN 

THE 

UNDISTORTED 

CRYSTAL 

R 1 (A°) 

e i 

*1 

R 2 (a°) 

e 2 

*2 

2.246 

67°05* 

-33°00 ' 

2.293 

142° 

03’ -36°59 ' 

(Rl 9 X 

(j^) and (R 2 © 2 ^2^ are 

the position coordinates of 


the first and the second nearest neighbours. 


* Reference [l] 

** Reference [4] 
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a 



(3 Fluorine ions (nearest neighbour) 

Fluorine ions (next nearest neighbour) 

• Y-ions 

+ Ions above the horizontal plane 
— Ions below the horizontal plane 

FIG. 3. 2 ; PROJECTED POSITIONS OF FLUORINE NEAR NEIGHBOURS OF 
A Y ION IN LiYF 4 SINGLE CRYSTAL. 
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TABLE 3.3 
IONIC RADII 


S.No. Element 


Ionic radius(A°) 


1 . 

Li + 

0.68 

2. 

y3+ 

0.893 

3. 

F“ 

1.33 

4. 

Eu 3+ 

0.950 


4 
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In the absence of the fluorine ions the symmetry of 
the unit cell would be D 0d . The two vertical reflection 
planes intersecting along the c-axis and the 2 two-fold 
rotation axes perpendicular to the c-axis can be easily 
visualized from Fig 3.1. The fluorine ions, however, do 
not possess these symmetry elements and the site symmetry 
of Y ions is lowered to with a 90° rotation about the 
c-axis followed by a reflection in a plane perpendicular 
to the c-axis. 

2.1 RARE EARTH DOPING IN LiYF ^ CRYSTALS 

Like other scheelites, LiYF 4 can be readily doped with 
rare earth ions. Since these crystals have compact struc- 
ture, interstitial substitution is very much unlikely. 

3+ + 

Therefore rare earth ions can replace either Y or Li 
ions. If the rare earth ions were to replace Li + , the 
resulting optical spectra would be very nearly unpolarized 
due to approximately cubic environment of the Li ions. The 
fact that the observed spectra are highly polarized suggests 
that the rare earths do not substitute for Li + ions. More- 
over, large differences in the ionic radii and ionic charges 
of rare earth and Li ions would make this substitution less 
likely. Conditions are much more favourable for a rare 
earth ion occupying a Y 3+ site on both these counts. These 
considerations show that LiYF 4 is a superior host for rare 
earth doping. 


I ' . V . • > 

a Jn«4vs 


AT* 
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3. EXPERIMENTAL SET-UP 

3 . 1 EXPERIMENTAL SET-UP FOR ABSORPTION MEASUREMENTS 

Figure 3.3 gives a schematic diagram of the set-up used 
for absorption measurements. Light from a white light 
source ’S' is focussed by the lenses L-^ and L ^ onto the 
sample mounted in a cryostat. The transmitted light is 
collected and focussed on the entrance slit of a double 
grating monochromator (GDM-1000) . The spectrally dispersed 
light from the GDM-IOOO is detected by the PMT (M12FC51) 
attached to the exit slit of the monochromator. The PMT 
signal is converted into a voltage signal by an I A - V A 
converter and amplified by a lock-in-amplifier (HR-8, PARC). 
The source 'S' is a 1000 watts tungsten-halogen lamp with 
linear filament (Osram). The colour temperature of these 
lamps is 2800K. The lamp is water cooled (Fig 3.4). For 
measurements at liquid nitrogen temperature a glass dewar 
was fabricated using cold finger for sample cooling. The 
cryostat and the sample holder are shown in Figure 3.5. The 
sample holder was made from a copper block and brazed to the 
cold finger. The dewar is fitted with 4 windows of Corning 
7059 glass. 

The spectral range of the monochromator from 7500 cm 1 
to 35000 cm""* is scanned in two orders. The instrument 
resolution is 0.5 cm -1 . A 25 Hz chopper and a photodiode 
are built in the monochromator for lock-in-detection. The 
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FIG. 3. 3 : SET-UP FOR ABSORPTION MEASUREMENT 





FIG. 3. 5 : THE CRYOSTAT AND THE SAMPLE HOLDER 
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wave number read out of the monochromator was calibrated 
using low-pressure Hg and Ne discharges, and Ar + ion laser 
lines. Additional details of the-GDM-1000 and its calibra- 
tion are given in appendix II. 

3 . 2 EXPERIMENTAL SET-UP FOR RECORDING FLUORESCENCE AND 
EXCITATION SPECTRA 

Fig 3.6 shows schematically the arranoement used to 
record fluorescence and excitation spectra. In this setup, 
a laser beam from a Coherent ring dye lase^ using R6-G dye 
or a Coherent (lnnova-100) 15-V7 Ar-ion laser, was incident 
on the sample and the emitted fluorescence was collected in 
the perpendicular geometry. The detection and the recording 
systems are identical to the absorption measurements. 

To select a and m polarizations a sheet polarizer was 
used. In order to compensate for the polarization characteris- 
tics of the monochromator, the spectral profile of the 
tungsten-halogen lamp was recorded by Inserting the sheet 
polarizer between the lamp and the entrance slit of the 
monochromator. The lamp profiles were recorded with the sheet 
kept in a and r positions and the necessary correction 
factors (I /I ) were obtained in the reaion of our interest. 

71 C5 

p^g 3.7 shows the o— polarized, it— polarized and unpolarized 
lamp profiles, and Fig 3.8 shows the dependence of the 
correction factors (1^/^) on wave number. The experimental 
data for these plots are listed in appendix II. 



Power 
Supply H.V. 
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FIG. 3. 6 : SET-UP FOR FLUORESCENCE MEASUREMENTS. 
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FIG. 3. 7 : RESPONSE OF •GOM-IOCO' TO POLARIZED LIGHT. 




1 (I-order) 



PLOT OF THE POLARIZATION CORRECTION FACTOR (T^/l 
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4 * EXPERIMENTAL DATA 

In this section, we present experimental data as 
obtained from absorption, fluorescence and excitation 
spectra. The wave numbers corresponding to the peak posi- 
tions of spectral lines were measured from the chart 
recording with reference to the positions of marking signals 
and corrected for monochromator calibration. The accuracy 
of measurement of the peak positions is within +1 cm""^ for 
sharp lines and within +2 to +5 for broad and weak lines. 

The wave number assignments were found to be always repro- 
ducible and were within this precision. The experimental 
parameters and conditions such as slit-width, and the 
sample orientation with respect to the incident radiation 
and the direction of observation were identical for o and it 
polarized spectra for any given group of transitions. 

Though we have tried our best to keep these conditions the 
same for recording the axial spectra as well, some departures 
are expected as the sample had to be physically rotated and 
because the axial spectra were usually recorded on different 
days. The sample temper ature was slightly higher than the liquid 
nitrogen temperature because it was not in direct contact 
with the liquid nitrogen. This may affect the results 
because spectral line positions in the LiYF 4 :Eu^ + are tempera- 
ture dependent. In fact, absorption spec tral lines at room 
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temperature are shifted upto 7 cm 1 from those at liquid 
nitrogen temperature. 

4.1 THE ABSORPTION DATA 

3+ 

The transverse and axial absorption spectra of tlYP^sEu 

were recorded at liquid nitrogen temperature in the spectral 

region 18500 cm 1 to 25500 cm” 1 . The absorption spectra are 

reproduced in Fig 3.9 to Fig 3.12. We could not see any 

absorption corresponding to 7 F - 5 D and 7 F - 5 D_ transitions 

oo o 3 

in the regions 17200-17300 cm” 1 and 24000 cm” 1 -24500 cm” 1 

respectively. A slit width of 70p and a scanning speed of 
"“1 

20 cm” /min were used for these measurements. Measured 
wave-number position of lines after correcting for the 
mono-chroma tor error are tabulated in Table 3.4. 

4.2 FLUORESCENCE DATA 

Information on the low lying energy levels can not be 
obtained from the absorption spectra. Fluorescence spectra 
maps the energy levels below the excited levels. Fluor- 
escence spectra of the sample were recorded at liquid 
nitrogen temperature in the transverse and axial sample 
orientations. For recording the transverse spectra, the 
sample was excited along the c-axis and the observations 
were made along one of the a-axes, while for the axial 
spectra, excitation was done along an a-axis and observa- 
tions were made along the c-axis. The slit width ranging 
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from 50 to 200|a was used depending upon the strength of the 

fluorescence of various groups of transitions. Commonly 

used scan speed was 4 cm ^/min in 2 nc * order and 2 cm”^/min 
s 

in the 1 order. The sample fluoresces when irradiated 

with any of the Ar + laser lines, but the fluorescence is 

strongest for the excitation at 457.9 nm. We have recorded 

the fluorescence with the 457.9 nm and with the 514.5 nm 

lines of the Ar + ion laser. The 457.9 nm excitation gives 

fluorescence from D^t Dj, and 3 D q excited levels while the 

514.5 nm produces fluorescence from the D, and D levels 

1 o 

5 

only. To excite the D Q level exclusively we used the R6-G 
dye laser lasing at 16938 cm’"'*' and 17630 cm - '*'. The resonance 

5 

excitation of D q at 17271 cm is extremely weak as the 

^F q -^D 0 direct transition is highly forbidden. To excita- 

tion of D level by the 17640 cm line of the dye laser is 
o 

phonon assisted while excitation by the 16938 cm”'*' line 

7 5 

corresponds to the F^- D q transition indicating enough popula- 

7 

tion of the F^ level even at liquid nitrogen temperature. The 
fluorescence spectra are reproduced in Fig 3.13 through 
Fig 3.28. Tables 3.5, 3.6 and 3.7 list the observed line 
positions for 457.9 nm, 514.5 nm and R6-G excitations 
respectively. The lines marked 'r* and v are Raman and 
vibronic transitions respectively. The lines with marks 'a 1 , 

'b' and 'c' are unidentified, Hg-lines and Ar + - lines 


respectively 
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4.3 EXCITATION DATA 

We have not made a detailed study of the excitation 

spectra. Whatever measurements were made in this respect 

helped us to identify the interfering vibronic lines. They 

also helped us in ascertaining the position of the D q level 
3+ 

of Eu indirectly. To obtain the excitation spectra we 

5 7 —1 

monitored the D - F. fluorescence at 16840 cm and scanned 

o i 

the R6-G dye laser. The peak intensity of the fluorescence 
signal was plotted with respect to the birefringent micro- 
meter reading which was then calibrated in wave numbers. This 
calibration is rather inaccurate with a maximum error of 
+6 cm - '*'. The excitation spectra are reproduced in Fig 3.29 
to Fig 3.32. Table 3.8 gives the position of the excitation 
peaks with the reported phonon— modes of LiYF^ pure crystal 
[6]. Fig 3.33 reproduces the micrometer calibration curve. 
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(The cm"^ scale is without calibration correction) 
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SPECTRA OF LiYF 4 :Eu 3+ AT LIQUID N 2 TEMPERATURE 
(The cm - '*' scale is without calibration correction) 
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TABLE 3.4 


ABSORPTION DATA FOR Eu 3+ IN LiYF ^ AT LIQUID NITROGEN TEMPERATURE 


ADOunriiun unm i l- u xin J-x i i ^ ni Ll^UlU IM1 IrtUbCN ItMFtnAiUnr: 

SNo. 

Line 

Posi- 

tion 

(cm -1 ) 

FWHM 

( cm -1 

Re la tive 
Intensity 
) (arbitrary 
units ) 

Transi- 
tion 
Mech . 

- S 4 -Sym- 
me try 
assign- 
ment 

Transition 





X I 

71 

ax 




1. 

25416 

34 

189 

141 

284 

TtE 

oE 

n-i5<2) 

ri-R< 3 > 

7p 0 - \ 

2. 

25258 a 

12 

17 

0 

16 

- 

- 

- 

3. 

25241 

14 

29 

0 

50 

oE 

r r r^( 2 ) 

?F o - 5l 6 

4. 

25027 

4 

0 

17 

0 

teE 

q-r 2 ( i) 

?F o - % 

5 

24953 

3 

5.0 

0 

8.0 

oE 

R-r 3 u) 

?P o - \ 

6. 

21543 

3 

0 

13 

0 

uE 

q-r 2 (2) 

?F o - \ 

7. 

21520 

3 

14 

0 

16 

crE 

r i" r 3 

?F o - 5 °2 

8. 

2l518 ?,a 

- 

- 

- 

- 

- 

- 


9. 

21450 

3 

0 

4 

0 

TtE 

p r r 2 d) 

7p o " 5d 1 

10. 

19043 

2 

6.0 

0 

0 

crM 

R-ri 

7 S 

F o - D 1 

11. 

190 42 s * a 

- 

- 

- 

- 

- 

- 

- 

12. 

19021 

2 

0 

7.0 

7.0 

TtM 

tI-«3 

7 f - 5 D. 
o 1 

13. 

19019 s ,a 

- 

- 

- 

- 

- 

«• 


14. 

186 94 a 

B 

0 

3 

- 

- 



15. 

18684 

4 

0 

4 

3 

itE 


7 5 

F 1 - D 1 
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a - refers to unidentified line 
s - refers to weak shoulder 
B - refers to broad and weak line 

n and cr indicate whether the electric vector of the 
radiation is | | or 1 "to the c crystal axis respectively. 

E and M indicate electric dipole and magnetic dipole 
transitions . 


-x- 





^ *n n»^u 
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20298 203?8 20458 20538 v 20618 206! 

cm" ^ 


FIG. 3. 14 : 5 D 2 - 7 F 2 FLUORESCENCE OF LiYF 4 :Eu 3+ AT LIQUID N 2 

TEMPERATURE 

(a and c indicate unidentified and Ar + lines 
respectively; The cm”^ scale is without calibratioi 
correction) 



Relative Intensity (I) 






FIG. 3. 16 ; 


j D 2 - F 4 FLUORESCENCE OF LiYF 4 :Eu AT LIQUID N 2 
TEMPERATURE 

(a and b indicate unidentified and Hg-lines 
respectively; The cm”^ scale is without calibration 
correction) 





Relative Intensity, (I) 
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TEMPERATIRE 

(a and b indicate unidentified and Hg-lines 
respectively; The cm” 1 scale is without calibration 
correction) 


f 






Keiative intensity 


axial -pol 



tt - pol 


cr-pol 


16488.5 


16588.5 

cm 1 


16688.5 


16788.5 


FIG. 3. 18 : V 7 F 6 flu0RESCENCE 0F LiYF 4 :E,j3+ AT LIQUID N, 
temperature 

(This group has not been analyzed? The cm - ^ 
is without calibration correction) 


scale 




Relative Intensity^! ) 



-1 

cm 

FIG. 3. 19 ; - 7 F AND 5 D.- 7 F. FLUORESCENCE OF LiYF . : E u 3 

loll ** 

AT LIQUID N 2 TEMPERATURE 

(The cm ^ scale is without calibration correc 





Relative Intensity, (I ) 



16968.5 17008.5 17068.5 17128.5 17181 

cm* 

FIG. 3. 21 : 5 D 1 - 7 F 3 FLUORESCENCE OF LiYF 4 :Eu 3+ AT LIQUID N 2 


TEMPERATIRE 

(The cnT 1 scale is without, calibration correction) 


Relative Intensity, (1 ) 



TEMPERATURE 

(The cm”^ scale is without calibration correction) 





cm 1 

FIG. 3. 23 : 5 D r 7 F 5 FLUORESCENCE OF LiYF 4 :Eu 3+ AT LIQUID 
TEMPERATURE 

(a indicates unidentified lines; The cm” 1 sc 
without calibration correction) 










Relative Intensity, (I ) 



152 6 8.5 15328.5 15388.5 15448.5 


cm 1 

FIG. 3.26 : 5 D o - 7 F 3 FLUORESCENCE OF LiYF 4 :Eu 3+ AT LIQUID N 2 

TEMPERATURE 

(The cm ^ scale is without calibration correction) 



Relative Intensity, (I ) 



FIG. 3. 27 : 5 D o ~ 7 F 4 FLUORESCENCE OF LiYF^Eu 34- AT LIQUID N 2 

TEMPERATURE 

(The cm ^ scale is without calibration correction) 



Relative Intensity, ( I ) 


1X1 



13208.5 13248.5 13288-5 v 13408.5 13448.5 13488.5 

cm " 1 

FIG. 3. 28 : 5 D - 7 F- FLUORESCENCE OF LiYF /1 :Eu 3+ AT LIQUID N 0 
O D 4 Z 

TEMPERATURE 

(The cm - '*' scale is without calibration correction) 
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TABLE 3.5 


FLUORESCENCE DATA FOR THE 457.9 nm EXCITATION OF Eu^tLiYF ^ AT 
LIQUID NITROGEN TEMPERATURE 


SNo. 

Obs . 

Line 

Posi- 

tion 

(cm -1 ) 

FWHM 

(cm -1 ) 

Relative 
Intensi ty 
(arbitrary 
units ) 


Transi- 

tion 

Mecha- 

nism 

S .-Sym- 
metry 
assign- 
ment 

Transi- 

tion 


I 

a 

F ax 

1. 

21570“’ r 

- 

- 

- 

- 

- 

- 

- 

2. 

21543. 

2.0 

9.2 

1.9 

0.4 

tcE 

r 2 <2)-n 

5 d - 7 f 

U 2 h o 

3. 

21520. 

2.0 

3 .6 

18.0 

10.5 

oE 

r 3 -tr 

5 n 7 p 
° 2 ~ F o 

4. 

21509? 

4.0 

2.4 

1.9 

2.0 

- 

- 

- 

5. 

21470?’ C 

2.0 

3.2 

2.9 

0.8 


- 

- 

6 . 

21463? 

2.0 

1.2 

0.0 

0.0 

- 

- 

- 

7. 

21458? 

2.0 

0.8 

1.4 

1.4 

- 

- 

- 

8 . 

21450. 

2.0 

10.0 

2.4 

0.8 

nE 

r^iui 

5 d 0 - 7 f 

2 o 

9. 

21408? 

3.5 

0.4 

3.4 

1.2 

- 

- 

- 

10. 

21209. 

3.0 

4.8 

1.2 

5.2 

■jtM 

r 2 ( 2 )-q 

5 7 

V F 1 

11 . 

21187. 

3.0 

6.0 

4.2 

2.4 

tcE , oM 

r 3 -u 

— 7 F 
U 2 1 

12 . 

21158? 

2.0 

3.6 

3.0 

1.8 

- 

- 

- 

13. 

21145. 

2.5 

7.0 

11.5 

19.4 

oE ,itM 


— 7 F 
u 2 ^1 

14. 

21117. 

2.5 

7.0 

7.2 

12.6 

icM , crE 

r 2 (D-E 

5 d - 7 f 
U 2 ^1 
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15. 

21090. 

2.5 

3.5 

0.4 

5.5 

itM 

r 3 -ir 

D 2 - F 1 

16. 

21048*? 

— 

0.0 

0.2 

0.0 

cM 

rail 

“*D — 7 F 
°2 ^1 

17. 

21021 . 

3.0 

3.0 

0.0 

0.0 

nE 


“*D - 7 F 
°2 h l 

18. 

20866 d . 

3.0 

48.6 

10.9 

15.2 

nE , oM 

- 

— 

19. 

20788*? 

4.0 

4.0 

8.3 

16.6 

TcM,aE 

— 

mm 

20 . 

20706? 

7.0 

4.9 

0.0 

0.0 

- 

— 


21 . 

20671*? * a 

4.0 

2.0 

0.0 

0.0 

- 

— 


22 . 

20629. 

2.5 

14.0 

16.3 

19.0 

nM , aE 

r 3 -qti) 

— ^F 

°2 b 2 

23. 

20589. 

3.0 

9.6 

1.6 

0.0 

tcE 

n-^ (i) 


24. 

20578*?’ 3 

4.5 

8.1 

0.0 

0.0 

- 



25. 

20567. 

4.0 

4.0 

5.2 

14.4 

itM, aE 

r 2 ( 2 )-r 3 


26. 

20545. 

4.0 

56.8 

11.7 

10.9 

nE , ( crM) E- [ 3 * 

— ”^F 

u 2 r 2 

27. 

20502 

4.0 

8.8 

21.8 

52.1 

itM, aE 

n - 13 

5 D - 7 F 

U 2 2 

28. 

20493? 

2.0 

20.4 

19.2 

5.9 

- 

. 


29. 

20474 . 

4.0 

12.8 

30.2 

72.3 

itM, crE 

r z w-r 3 

u 2 r 2 

30. 

20455 . 9 


0.4 

0.4 

0.0 

- 



31. 

20390*? 

. 

0.3 

0.3 

0.0 

aM 

r^(2)-fj(2) V 1 

32. 

20369*? 

B 

1.0 

0.8 

0.9 

cjE,tcM 

r 2 ( 2 )-ri 

and 

r 3 -i2(2) 

5 n -'f 
° 2 h 2 

33. 

20344°? 


0.5 

0.8 

1.5 

aE , (tcM) r 3 - rj 

_"^f 

u 2 2 
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34. 

20271. 

3.0 

4.8 

4.2 

0.0 

nE 

rj(D-ri 

5 d - 7 f 
u 2 t 2 

35. 

19669. 

3.0 

4.2 

3.2 

9.2 

nM , cE 

r 2 ( 2 )-i;(i) 5 D 2 - 7 r 

36. 

19659. 

2.5 

7.5 

3.2 

12.5 

7tM,crE 

fs-Jad) 

5 d - 7 f 

2 ^3 

37. 

19646. 

3.0 

2.4 

3.8 

2.3 

nE f aM 


V ?F 3 

38. 

19617. 

2.5 

69.0 

11.3 

21.1 

nE 

and 
crE t nM 

q-rjd) 

and 

r 3 - p 

V ?F 3 

39. 

19605. 

2.5 

10.5 

9.0 

22.6 

TtM, ctE 


5 d - 7 f 

2 r 3 

40. 

19590. 

3.0 

7.8 

12.6 

7.5 

hM, <tE 
and 
a . M 

Y 2 (2)-V 3 (2) 5 D 2 - 7 F. 
and 

r 2 u)- gd) 

41. 

19576 . 

3.0 

37.8 

20.7 

34.1 

TtM, cE 
and 
crM 

r^(i)-^(i) _ 7 

r r if d 

42. 

19566. 

2.5 

18.5 

7.5 

4.3 

irE , crM 

<> 


43. 

19548. 

2.5 

49.0 

3.6 

5.0 

nE 

Y 2 U)-Y 1 

5 d - 7 f 

u 2 t 3 

44. 

19525. 

2.5 

8.0 

1.8 

10.9 

TtM, (ctE)^- V 3 (2) 

5 n _ 7 f 

U 2 b 3 

45. 

19504? 

- 

0.0 

0.5 

0.3 

oM 

P 2 (2)-[J(2) 5 D 2 - 7 F. 

46. 

19496 . 

3.0 

10.2 

8.3 

19.6 

tcM , ctE 

r 2 (2)-i^(2) 5 d 2 - 7 f s 

47. 

19480. 

7.0 

7.0 

1.4 

8.1 

nM 

r 3 -rj(2) 

R 7 

— / F 

U 2 h 3 

48. 

19439? 

6.0 

2.4 

0.0 

0.0 

nE 

ri-r 2 <2) 

— 7 F 

u 2 r 3 

49. 

19042. 

2.0 

6.0 

45.6 

10.0 

aM 

1-1? 

5d 1- 7f o 

50. 

19020. 

2.5 

96.0 

15.0 

35.0 

tcM,crE 

*3-17 

5d 1- ?F o 

51. 

18912? 

- 

- 

- 

- 

— 

*3-*I (1 > 

V ?F 4 
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52. 

18842“ 

- 

- 

- 

- 

- 

ryij-fjd) 5 d 2 - 7 f 

53. 

18793“’ a 

- 

- 

- 

- 

- 

— — 

54 . 

18709. 

2.0 

11.6 

25.2 

24.6 

nM , crE 


55. 

18687. 

2.0 

36.4 

5.4 

2.5 

nE , csM 

15 - 13 5d i- ?f i 

56 . 

18666. 

6.0 

12.0 

2.7 

12.7 

rtM, aE 

V ?F 4 

57. 

18636. 

5.0 

0.0 

4.5 

14.4 

aE 

r 2 (i)-r 3 (i) 5 d 2 - 7 f, 

58. 

59. 

18619? 

18612?’ 03 

- 

- 

— 

— 


5 d 2 - 7 f 4 

n- n v 7 5 

60. 

18602? ,a 

- 

- 

- 

- 

- 

- - 

61. 

18592. 

3.0 

15.0 

33.5 

22.2 

nM, aE 

11 - tr 5 °i- 7f i 

62. 

18580?’ a 

- 

- 

- 

- 

- 

- 

63. 

18574?’“ 

- 

T* 

- 

- 

- 

n-i^u) 5 d 2 _ 7 f 4 

64. 

18565“ 

- 

- 

- 

- 

- 

£(Z)-£(2) 5 D 2 - 7 F 4 

65. 

18542. 

7.0 

11.2 

31.9 

62.1 

aE , 7i M 
oM 

r 3 -q<2) 

r^fD-gd) 

66. 

18528. 

3.5 

32.2 

7.6 

35.6 

tcM, aE 

fj(2)-ij(2) 5 D 2 - 7 F 4 

67. 

18505 

3.5 

3.5 

5.6 

12.8 

aE,7CM 

r 3 -d(2) 5 d 2 - 7 f 4 

68. 

18492? 

3.5 

22.4 

1.2 

0.0 

- 

- 

69. 

18481? 

4.0 

8.0 

12.6 

18.6 

- 

- 

70. 

18471?’“ 

- 

-r 

- 

- 

- 

5< ll -G<2> V ?F 4 

71. 

18467. 

B 

1.8 

3.6 

6.0 


fl-fj(2) 5 D 2 - ?F 4 
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72. 

18455?’ a 

73. 

18436. 

74. 

18425?’ a 

75. 

18407?’ a 

76. 

18394? 

77. 

18379?’ a 

78. 

18373? 

79. 

18358? 

80. 

18332? 

81. 

18303? 

82. 

18247? 

83. 

18217? 


84. 18182?’ a 

85. 18179?’ a 

86. 18167?’ a 


CD 

• 

18153? 

- 

CO 

CD 

18144? 

5.5 

89. 

18131. 

3.0 

90. 

18067, 

3.5 

91. 

18043. 

3.5 


0.6 

0.4 

0.0 

- 

- 

- 


• 21.0 

4.8 

11.5 

aE , 

ttM 

R-(l(3) . 7 




oE , 

n'A 

qu)- 

-R(3) °2~ 

0.0 

0.8 

1.7 

— 



O 

1.6 

1.3 

0.0 





2.6 

1.7 

3.8 

- 


WO) 

0.0 

1.3 

3.0 

- 




1.4 

0.0 

3.0 

_ 




1.0 

0.0 

2.2 





1.4 

18.5 

3.1 





0.8 

0.4 

1.6 

— 




1.4 

0.0 

0.0 

— 




0.0 

13.8 

2.9 

- 


- 

- 

1.4 

0.0 

0.0 




— 

2.8 

0.0 

0.0 

- 


- 


0.2 

0.0 

0.0 

7lE 

ri- 

-gu) 

5 d i - 7 f 2 

38 . 5 

0.0 

0.0 

- 


- 


50.4 17.3 54.3 

, aE 

G- 

•ijd) 

5d i- ? f 2 

30.1 14.0 43.0 Tty 

f aE 

n- 

G 

5d i- 7 f 2 

12.9 37.5 i 

3.2 crM, 


r 3 - 

G 

5 d,- 7 f„ 


F 
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92. 

17891? 

— 

4.0 

0.0 

0.0 

TtE 

II-|J( 2 ) 5 D r 7 F 2 

93. 

17870. 

7*5 

142.5 

78.0 

156.0 

ttM, ctE 

rj-ij( 2 ) 5 Di - 7 f 2 

94. 

17845! 

- 

6.0 

3.6 

7.9 

nM, oE 

7-7 5 d x - 7 f 2 

95. 

17813?’ a 

- 

1.6 

1.1 

1.6 

- 

- - 

96. 

1777l! 

- 

12.4 

8.4 

15.6 

- 

- - 

97. 

17745? 

- 

8.2 

6.2 

10.8 

■rcM,crE 

7(2)- 7<i) 5 d 2 - 7 f, 

98. 

17719? 

- 

9.4 

3.4 

5.7 

7tM,crE 

^-1^(1) 5 D 2 - 7f 5 

99. 

17691? 

- 

7.4 

0.0 

4.5 

ttE 

H-ija) 

100. 

17678? 

- 

4.0 

2.9 

5.9 

itM, 0 E 

5 d 2 - 7 f 5 

101. 

17605. 

4.0 

20.8 

1.2 

1.0 

TtE 

rj(D- rjd) %- 7 f, 

102. 

17533. 

9*0 

9.0 

0.0 

0.0 

tcE 

t2(2D-fr< 2 ) X- 7 ^ 

103. 

17480? 

- 

0.0 

4.0 

6.4 

crE 

2 ) 5 d 2 - 7 f 5 

104. 

17451. 

B 

0.0 

8.2 

11.9 

a E 

rj(i)-rj(2) 5 d 2 - 7 f c 

105. 

17443. 

B 

20.5 

0.0 

0.0 

tiE 

v 2 M-rj(2) 5 d 2 - 7 f 5 

106. 

17428. 

- 

0.0 

1.6 

2.8 

cE 

r r ij(3) 5 d 2 - 7 f 5 

107. 

17400? 

- 

0.0 

0.8 

1.4 

crE 

q(D-go) 5 d 2 - 7f s 

108. 

17376? 

- 

0.0 

0.6 

0.6 

0 

nci>-r^C2) V 7 f 5 

109. 

17335. b 

B 

3.9 

0.8 

1.4 

- 

- - 

110. 

17272** 

B 

9.0 

4.8 

5.0 

- 

0- 7 X-X 

111. 

17253? 

B 

7.7 

6 .2 

5.8 

mm 
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16714 


133. 16686. 

134. 16667. 

135. 16656. 

136. 16642. 

137. 16640. 

138. 16625. 

139. 16607. 

140. 16600. 

141. 16583. 

142. 16550. 

143. 16511. 

144. 16508. 

145 16490. 

146. 16486. 

147. 16470t 



94 


148. 

16436“ 

- 

0.0 

CM 

• 

0.0 

c jW 

n-qu) 

5 7 

q- ^ 

149. 

16413. 

5.0 

57.4 

87.0 

154.4 

irM , oE 

q-iid) 

5d 1- 7f 4 

150. 

16382. 

3.0 

122.4 

17.5 

21.5 

tiE 

Prfjd) 

5d o- 7f 2 

151. 

16295. 

4.5 

71.8 

75.6 

169.7 

txM, aE 

q-q 

- 7 F 
o r 2 

152. 

16229. 

6.0 

38.9 

34.0 

82.6 

7i M , dE 

q-qaj 

5 d - 7 f 

1 ^4 

153. 

16208. 

4.5 

160.6 

28.0 

42.6 

tiE , crM 

q~qd) 

5 d - 7 f 

U 1 ^4 

154. 

16149. 

6.0 

42.0 

25.2 

71.2 

oE , tcM 

q-q^) 

— 7 F 

U 1 h 4 

155. 

16138. 

5.0 

283.8 

38.0 

73.9 

tiE , crM 

n-rjd) 

5 7 

D - F 
U 1 *4 

156. 

16121. 

14.0 

- 

- 

- 

tcE 

q-r?d) 

5 °o- 7f 2 

157. 

16116. 

B 

- 

- 

- 

aE 

ri-Rd) 

V 7 % 

158. 

16089“ 

- 

- 

- 

- 

- 

n-ir 

\- ?F 2 

159. 

16064“ 

- 

5.8 

0.0 

0.0 

itE 

q-gts) 

5 d - 7 f 

U 1 r 4 

160. 

16042. 

10.0 

191.4 

184.0 

331.9 

nM, crE 


- 7 F 

U 1 M 

161. 

16032! 

B 

- 

- 

- 

itM, crE 

ri-i3<2) 

- 7 F 

U 1 ^4 

162. 

16009. 

6.0 

236.6 

34.8 

42.8 

itE , crM 

r 3 -i3<2) 

- 7 F 

U 1 t 4 

163. 

15940. 

17.0 

153.0 

129.2 

289.4 

itM, aE 

r3-ir< 3 > 

\-\ 

164. 

15411. 

3.0 

70.2 

7.8 

10.1 

itE 

fl-fld) 

\-\ 

165. 

15398. 

3.0 

16.6 

21.0 

45.0 

crE ,7tM 

q-qw 


166. 

15368. 

3.0 

0.0 

12.6 

0.9 

a M 

q-q 

O o 

167. 

15317. 

3.3 

29.7 

40.3 

81.2 

aE,itM 


5 ° 0 - 7f 3 

O o 
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168. 

15283? 

10.0 

0.0 

26.0 

37.9 

— 

- 

wm- 

169. 

15245. 

4.0 

15.8 

36.8 

72.1 

crE ,-rcM 

ri-rju) 

U 1 r 5 

170. 

15236.. 

5.5 

43.1 

58.3 

115.0 

0 E , nM 

17 - 17(2) 

5d o- 7 F 3 








13 - 1? 1 ’ 

5 d - 7 f 

1 5 

171. 

15223. 

7.0 

68.6 

22.4 

3.8 

nE, (crM) j^-fjCl) 

5 n 7 p 
D l' F 5 

172. 

15214“ 

- 

0.0 

3.0 

4.4 

- 

- 


173. 

15203“ 

- 

0.0 

0.0 

3.5 

- 

- 

- 

174. 

15191? 

10.0 

38.0 

46.0 

65.1 

- 

- 

- 

175. 

15177. 

6 .0 

11.3 

25.2 

54.7 

ctE ,-icM 

r 3 -ri(i) 

“*D - 7 F 

U 1 D 

176. 

15044“ 


0.0 

0.4 

0.6 

crE 

n-rjts) 

“*0 - 7 f 

U 1 0 

177. 

15022. 

5.0 

49.2 

10.0 

0.0 

nE , crM 

r 3 -r 3 (2) 

- 7 F 

U 1 r 5 

178. 

15014. 

11.0 

0.0 

26.4 

47.9 

aE 

r 3 -rc<2> 

- 7 F 

U 1 r 5 

179. 

14994. 

9.0 

14.4 

28.8 

56.4 

crE ,7tM 

n-r 3 ( 3 ) 

- 7 F 

U 1 r 5 

180. 

14971. 

7.0 

145.6 

30.6 

41.0 

tcE , crM 

13-13(3) 

5 7 

D l- F 5 

181. 

14963“' a 

- 

0.0 

2.0 

3.1 

- 

- 

— 

182. 

14457 

4.5 

71.3 

74.0 

87.3 

icM, crE 

fl-r 3 u) 

5 d - 7 f. 

0 4 


183. 14400. 9.0 4.9 43.2 1.3 crM 5 D q - 7 F 4 

184. 14367 . 5.0 85.1 10.0 8.3 itE,aM ( 1 ) 5 D q - 7 F 4 

185. 14293. 9.5 526.7 70.2 61.4 irE ^-r3<2) \- 7 f 4 

107.5 87.6 122.9 nM.oE n-G< 2 > v 7f 4 


186 


14259 


7.0 
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* - Relative intensities have meaning within a aroup . 

a,b, c: - Refers to unidentified, Hq-line (tube light)and Ar + 
lines respectively 
d - Pr^ + lines 

r - Refers to Raman lines 

V - Refers to vibronic lines 

s - Refers to soulder 

w - Refers to very weak lines 

B - Refers to broad lines 

** - Assignment of 17272 cm - '*' is not certain 

+ - ++ This group of lines from 16816 cm~^ to 16470 cm ^ 

5 7 

may correspond to D^~ F^: We could not analysed 
these . 
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TABLE 3.6 

FLUORESCENCE DATA FOR THE 514.5 nm EXCITATION OF Eu 3+ : LiYF^ 
AT LIQUID NITROGEN TEMPERATURE 


SNo. Obs . 
Line 
Posi- 
tion 

FWHM 
( cm -1 ) 

Relative * 
Intensity 
(arbitrary 
units) 

Transi- 

sition 

Mecha- 

nism 

S .-Sym- 
metry 
assign- 
ment 

Transi 

tion 

(cm -1 ) 


i I 

n a 

Lx 




1. 

19043. 

2.0 

1.2 

5.7 

0.8 

ctM 

a-n 


2. 

19021. 

2.0 

18.0 

4.3 

13.0 

tcM, crE 

HrE 


3. 

18710. 

3.0 

17.4 

29.5 

41.7 

nM , aE 

u-u 

U 1 h l 

4. 

18687. 

2.5 

33.3 

5.1 

4.5 

tie, (om) r^-n^ 

— 7 F 
U 1 U 

5. 

18593. 

3.0 

11.4 

29.7 

46.3 

m f oE 

[Ml 

5 d - 7 f 

U 1 t l 

6. 

18494? 

2.5 

5.0 

2.0 

3.8 

- 

- 


7. 

18467? 

2.5 

3.0 

0.0 

0.0 

- 

- 

- 

8. 

18435? 

2.5 

1.0 

2.0 

2.9 

- 

- 

- 

9. 

18419? 

2.5 

1.0 

0.0 

0.0 

- 


- 

10. 

18382? 

2.5 

3.5 

0.0 

0.0 

- 

- 

- 

11. 

18153? 

- 

0.8 

0.0 

0.0 

irE 

q-qa) 

~*D — 7 F 

1 r 2 

12. 

18132. 

3.5 

44.8 

15.1 

46.8 

tcM, aE 

i) 

- 7 F 

U 1 r 2 

13. 

18067. 

4.0 

30.0 

11.8 

31.9 

itM, aE 

r r g 

- 7 F 
U 1 ^2 
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14. 

18044. 

4.0 

8.0 

37.7 

4.1 

itE , crM 

r 3 -r 3 

- 7 F 

U 1 t 2 

15. 

17892. 

4.0 

1.6 

0.0 

0.0 

tiE 

r r ri(2) 

5 d - 7 f 

L»1 t 2 

16. 

17869. 

4.5 

9.9 

2.6 

25.4 

nM,crE 

r 3 -!3fe> 

- 7 F 
U 1 t 2 

17. 

17842. 

9.0 

3.6 

0.0 

2.9 

1 JiM 

r 3 -fi 

- 7 F 

U 1 tr 2 

18. 

17185. 

3.5 

8.7 

0.1 

0.6 

tiE 

rr-nti) 

5 d - 7 f 

U 1 o 

19. 

17171. 

3.5 

17.4 

24.5 

38.6 

crE ,tcM 

rr-Ed) 

5 d - 7 f 

U 1 0 

• 

o 

CNJ 

17162. 

3.5 

29.9 

44.8 

73.1 

crE ,nM 

0-8d> 

5 d - 7 f 

U 1 ^3 

21. 

17149. 

4.0 

81.1 

8.0 

9.9 

tcE 

8-0 d) 

5 — 7 F 

\ F 3 

22. 

17119. 

3.5 

16.8 

20.3 

24.7 

crE ,-jiM 

8-fT 

5d 1- ?F 3 

23. 

17091. 

3.5 

40.6 

42.7 

95.5 

crE ,itM 

8-0< 2 > 

5 °1- 7f 3 

24. 

17069. 

3.5 

32.5 

7.0 

7.6 

nE 

r 3 - r 3 ( 2 ) 

5d 1- 7f 3 

25. 

17006 . 

6.0 

80.6 

6.6 

13.0 

itE 


— 7 F 

U 1 h 3 

26. 

16983. 

6.0 

10.0 

18.0 

29.3 

tcE,(itM) 

0-HC2) 

5 d - 7 f 

U 1 h 3 

27. 

16938. 

3.5 

241.9 

52.5 

160.0 

n:M,crE 

8-8 

5d o- 7f 1 

28. 

16932?’ a 

- 

- 

- 

- 

- 

- 

- 

29. 

16842. 

3.5 

25.0 

64.4 

14.0 

0 M 

8-8 

5 d - ?F i 

30. 

16411. 

5.5 

0.0 

5.5 

3.5 

crE 

O-Od) 

— 7 F 
U 1 *4 

31. 

16381. 

4.5 

25.2 

1.8 

4.0 

irE 

8-8 d) 

5d o- 7f 3 

32. 

16294. 

5.5 

48.8 

71.5 

112.5 

crE ,irM 

8-r 3 

V ?F 3 

33. 

16230*? 

- 

0.0 

.2 

.8 

nM f crE 

8-8 d) 

5d 1- ?F 4 
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34. 

16208. 

4.0 

14.0 

2.0 

5.5 

tiE , oM 

r 3 -r 3 (D 

5 d - 7 f 

U 1 ^4 

35. 

16138. 

4.0 

13.2 

0.2 

4.1 

itE , crM 

q-r^D 

5 d - 7 f 
U 1 M 

36. 

16149“ 

B 

0.0 

0.0 

0.0 

crE ,7iM 

r 3 -q(2) 

- 7 F 

U 1 M 

37. 

16121“ 

B 

1.9 

0.1 

0.0 

tcE 

r r £(2) 

5 d - 7 f 3 

00 

CO 

16116“ 

B 

1.9 

0.5 

1.8 

icM, aE 

p 3 -nu> 

“*D - 7 F 

U 1 ^4 

39. 

16042“ 

8.0 

10.0 

6.4 

27.8 

it M , crE 

ri-r^ca) 

5 7 

— F 

U 1 

40. 

16009. 

6.0 

73.1 

6.0 

17.0 

7i E , orM 


- 7 F 

U 1 t 4 

41. 

15939“ 

B 

0.0 

0.2 

4.8 

crE 

r 3 -r 1 o) 

5[) 1- 7f 4 

42. 

15411. 

6.0 

144.0 

24.0 

14.4 

irE 

r r p 2 d) 

V ?F 3 

43. 

15399. 

5.0 

36.8 

37.0 

39.0 

tiM, crE 

q-rj(D 

V ?F 3 

44. 

15369. 

4.5 

0.0 

20.7 

1.2 

eM 


V ?F 3 

45. 

15318. 

5.0 

40.0 

76.0 

72.7 

tiM, crE 

q~r5(2) 

5d o‘ 7f 3 

46. 

15270? 

5.0 

0.0 

2.0 

0.0 

- 

- 

- 

47. 

14458. 

6.0 

48.4 

67.2 

132.2 

crE ,nM 

r>r3(u 

\~\ 

48. 

14366. 

6.0 

82.9 

10.8 

12.1 

tiE 

rrfld) 

5 d 0- ?F 4 

49. 

14293 ; 

10.0 

528.0 

70.0 

81.8 

■hE 

r r r 2 ( 2) 

5 d - 7 f 4 

o 4 

50. 

14260. 

7.0 

61.4 

81.2 

175.5 

crE , tiM 

ri-r 3 (2> 

5 d - 7 f 4 

o ■ 4 


a - refers to unidentified lines 
to - refers to very weak lines 
s - refers to shoulder 
B - refers to broad lines 

* - Relative intensities have meaning within a group. 
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TA BLE 3.7 

FLUORESCENCE DATA FOR THE R6-G EXCITATION OF Eu 3+ ?LiYF 1 AT LIQUID 
NITROGEN TEMPERATURE 


SNo 

. Obs . 
Line 
Posi- 
tion 

(cm -1 ) 

FWHM 

(cm - '*' 

Relative 
\ Intensity* 
(arbitrary 
units) 


Transi- 

tion 

Mech. 

54-Sym- 

metry 

assign- 

ment 

Transi- 

tion 



I** 

ax 

1. 

17292? 

B 

0.0 

1.2 

0.0 

- 

- 

- 

2. 

17272. 

B 

1.0 

2.4 

0.0 

crM 

r\-q 

0 0 

3. 

17250? 

B 

0.0 

2.2 

0.0 

- 

- 

- 

4. 

17232? 

B 

0.0 

4.0 

0.0 

- 

- 

- 

5. 

16 938. 

2.3 

121.5 

45.0 

58.1 


q-q 

5d o- ?F 1 

6 . 

16933? 9 a 

2.3 

8.1 

1.8 

4.1 

- 

- 

- 

7. 

16841. 

2.3 

18.4 

28.8 

5.4 

crM 

q-q 

V 7f i 

8. 

16382. 

3.5 

22.0 

2.1 

1.5 

tiE 

q-q (i) 


9. 

16295. 

3.5 

43.3 

30.8 

64.4 

itM,crE 

q-q 

V ?F 2 

10. 

16122. 

9.0 

36.7 

5.4 

0.0 

tiE 

q-q (2) 

5d o- 7f 2 

11. 

15411. 

3.0 

54.0 

5.4 

1.8 

itE 

q-qu) 

— *^F 
u o t 3 

12. 

15398. 

3.0 

16.6 

11.4 

5.1 

cM,itE 

q-r 3 d) 

— ^F 
u o o 

13. 

15369. 

3.0 

0.0 

7.5 

0.0 

oM 

q-q 


14. 

15317. 

3.0 

18.2 

23.4 

8.7 

cM,nE 

q-q( 2) 

5d o- ?F 3 
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15. 

14664. 

6.0 

0.0 

6.0 

0.0 

a M 

q-qo.) 

V ?F 4 

16. 

14584? 

18.0 

0.7 

9.7 

0.9 

- 

- 

- 

17. 

14558? 

B 

0.0 

0.3 

0.0 

- 

- 

- 

« 

00 

rH 

14538? 

7.0 

0.8 

5.0 

0.3 

- 

- 

- 

19. 

14458. 

4.0 

14.1 

14.4 

18.8 

itM , aE 


V 7 f 4 

20 . 

14400. 

7.0 

0.8 

5.0 

0.3 

a M 

fl-qco 

5 d o- ?F 4 

21 . 

14366. 

4.0 

16.3 

0.2 

1.5 

tiE 

q-qu) 

5 °o- ?F 4 

22 . 

14294. 

8.0 

119.1 

8.0 

4.4 

nE 

ri-q< 2 > 

5 d - 7 f, 

0 4 

• 

CO 

CM 

14259. 

6.0 

11.5 

15.6 

23.2 

7cM,ctE 

n-q< 2 ) 

5 d - 7 f 4 

0 4 

24. 

13485. 

5.0 

55.5 

10.0 

0.0 

icE 

q-qa) 

5 [) o- 5f 5 

25. 

13472. 

4.0 

34.7 

62.4 

9.0 

oE 

q-qti) 

5 D 

O D 

26. 

13271. 

8.0 

20.8 

30.4 

0.0 

uE 

r r p 3 ( 2 ) 

X- 7f , 

27. 

13219. 

B 

0.0 

0.6 

0.0 

oE 

Hl“ [ 3 ( 2 ) 

5 ^o- 7f 5 

r - 

refers to 

Raman 

line 






a - 

refers to 

unidentified line 





s - 

refers to 

weak 

shoulder 






B - 

refers to 

broad 

and weak 

line 





* — 

Relative 

intensities have 

meaning 

within a group 



-**- The laser power used for axial spectrum was different from that 
for the other polarizations. 
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Micrometer Reading ( jj ) 


FIG. 3. 31 : EXCITATION SPECTRUM OF LiYF 4 :Eu 3+ AT LIQUID N 


TEMPERATURE 



Micrometer Reading(ju) 

FIG. 3. 32 .* EXCITATION SPECTRUM OF LiYF 4 :Eu 3+ AT LIQUID N 


TEMPERATURE 
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•e 

o 


FIG. 3. 33 : DYE LASER CALIBRATION CURVE 
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TABLE 3.8 


PEAK 

POSITIONS IN 

! THE EXCITATION SPECTRA OF Eu 3 " 1 " : LiYF* AT 

LIQUID NITROGEN TEMPERATURE AND 

ENERGY OF THE CORRESPONDING 

PHONONS 

SNo. 

Peak 

Positions 

(cm' 1 ) 

Difference 
from 5 d 0 - 7 F ( 
transi tion 

at 17270 cm"' 

Difference 

from 

L WtP 

transition 
at 16936 cm” 1 

Reported* 
Phonon . 

energy (cm” ' 

1. 

17718 

448 

782 

446Eg 

2. 

17702 

432 

766 

425Ag/427Eg/ 





424Eu 

3. 

17670 

400 

734 

396Au 

4. 

17662 

392 

726 

382Bg 

5. 

17624 

354 

688 

368Eg/367Eu 

6. 

17608 

338 

672 

329Bg/329Eg/ 





326Eu 

7. 

17568 

298 

632 

292/303Eu 

8. 

17552 

282 

616 

283AU 

9. 

17500 

230 

564 

- 

10. 

17492 

222 

556 

224Au 

11. 

17474 

204 

538 

199Eg 

12. 

17464 

194 

528 

195Au 

13. 

17420 

150 

484 

153Eg 
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14. 

17410 

140 

474 

143Eu 

15. 

17394 

124 

458 

- 

16. 

17378 

108 

442 

446 Eg 

17. 

17354 

86 

420 

425 Ag / 427Bg / 424Eu 

18. 

17334 

64 

398 

396Au 

19. 

17320 

50 

384 

382Bg 

20. 

17292 

22 

356 

368Eg/367Eu 

21. 

17276 

6 

340 

- 

22. 

17270 

0 

334 

- 

23. 

17250 

20 

314 

- 

24. 

17230 

40 

294 

292Eu 

25. 

17218 

52 

282 

282Bg/283Au 

26. 

17204 

66 

268 

264Ag 

27. 

17196 

74 

260 

264Ag/252Au 

CO 

CM 

17158 

112 

222 

224Au 

29. 

17144 

126 

208 

199E 3 

30. 

17124 

146 

188 

143Eu 

31. 

17104 

166 

168 

173Eu/l77Bg 

32. 

17060 

210 

124 

- 

33. 

17040 

230 

104 

224Eu 
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34. 

17020 

250 

82 

248/2 52 Bg/Au 

35. 

17004 

266 

66 

264Ag 

36. 

16 972 

298 

34 

292/30 3Eu 

37. 

16936 

334 

0 

. _ 


* - Reference [6]. 
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CHAPTER IV 


RESULTS AND DISCUSSIONS 

In this chapter we present the interpretation of the 
experimental data given in Chapter III. The spectra of the 
LiYF 4 :Eu^ + consist mainly of groups of sharp lines. These 
lines correspond to transitions among the Stark components 
belonging to various J-multiplets . The center of gravity of 
a certain group of lines gives the energy of the so called 
'free ion level'. The symmetrv assignments to Stark levels 
are made using E.D. and M.D. selection rules. The relative 
intensities of spectral lines in tc , o and axial soectra help 
decide the nature of the transition. Lines with comparable 
intensities appearing in o(n) and axial spectra are associated 
with ED(MD) transitions. Transitions are labelled as m-E.D., 
o-E.D., tc-M.D. and c-M.D. where o(ti) indicates that the 
electric vector of the emitted or absorbed radiation in the 
transition is perpendicular (parallel) to the C-axis of the 
crystal and E.D.(MD) implies that the transition is electric 
dipole (magnetic dipole) in nature. The following general 
rule is helpful in deciding the nature of observed transitions. 

A pure it -polarized E.D. transition must appear only in 
the it-polarized spectrum, a pure o-polari zed E.D. transition 
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must appear in a and axial spectra but not in n: -polar ization . 
While a pure n-polarized M.D. transition must aopear in tu 
and axial-spectra but not in 0 -polarization, and a pure 0 - 
polarized M.D. transition must aopear only in the 0 polar- 
ization. Mixed transitions will aopear in both (% and 0 ) 
the polarizations. Their relative intensities determine the 
dominant polarization and intensity in the axial-spectrum 
ditermines the nature of the transition. The criteria for 
mixed transitions are as follows. 

I_<I XI „ = 0 -E.D. (dominant), 11 . M.D. 

71 CJ 0 X 

I <^<1 = it -M.D. (dominant), 0 -E.D. 

I >1 >1 = 0-M.D. (dominant), tc-E.D. 

I >1 >1 v = tx-E.D. (dominant), 0 -M.D. 

71 O 0 X 

A comparison of E.D. and M.D. selection rules shows that the 
combinations like ( 0 .E.D. + 0 .M.D.) and (it. E.D. + it. M.D.) 
are not allowed. 

1. ABSORPTION SPECTRA 

The absorption spectrum directly gives the positions of 

3 + 

the excited energy levels. The ground state of Eu is a 
singlet "^F with J=0. The "^F, multiplets with J=l- 6 lie 
next in energy. A survey of literature on the spectra of 

Eu^" 1 " in various hosts and solutions reveals that the next 

b -1 

hiaher level D is expected to lie around 17270 cm 

^ o 

7 7 5. 

above F * but we could not observe the F - D transition 
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in- absorption. This transition is hiohly forbidden by 

several selection rules. The J=0 — > J' =0 transition is 

forbidden by E.D. and M.D. selection rules. F and D 

o o 

states both transform as representations of S^ point 
group and the ["* — transition is forbidden by the crystal 
field selection rules. We now describe the groups of levels 
observed in the absorption spectrum in the spectral range 
18600 cm -1 to 25500 cm" 1 . 

7 F i- 5d i absorption 

The absorption spectra (Fig 3.11 and 3.12) show two lines 

at 18684 and 18694 cm 1 . The line at 18684 cm 1 is weak 

— 1 

but sharp while the 18694 cm line is very weak and broad. 

The literature survey and analysis of the rest of our experi- 
mental data indicate that these lines do not correspond to the 

7 5 -1 

F q — transition. The line at 18684 cm ' present in 

n-oolar ization , absent in cr-polarization , and extremely weak in 
the axial corresponds to an E.D. transition and satisfies 
the polarization and energy difference requirements for the 

: X<rp transition. But we do not have any expla- 
nation for the very weak and broad line at 18694 cm" 1 . 

7 f - 5 d, absorption 

The J=1 multiplet splits into two Stark levels trans- 
forming as p x and (Table 2.3). For this qrouo, we have 
observed two sharp and strong lines at 19043 and 19021 cm" 1 
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and two very weak shoulders at 19042 and 19019 cm - '*' respect- 
ively (Figs 3.11 and 3 . 12 ). The line at 19043 cm - '*' corresponds 
to a c-polarized M.D. transition as it is strona in the 
cr-polar ization but absent in % and axial spectra. We can 

assign this line to the 7 F o (rp ( Fj ) transition. The 

line at 19021 cm -1 is strong in n and axial spectra but 
extremely weak in the o-polarization . This indicates that 
it is a n-polari zed magnetic dipole transition with a small 
admixture of o-E.D. character. We can attribute this line 

to the 7 F (r\) 5 D x (rp transition. The weak shoulders 

remain unexplained. 

7 f - 5 d^ absorption 

— o 2 . 

Theoretically we expect 4 lines for this group according 
to the decomposition D* ■ iq + 2r 2+ ir 3> There are three 

fairly strong and sharp lines at 21543 cm - ' 1 ', 21520 cm "* , 

21450 cm - '* and a weak shoulder at, 21518 cm - '*, associated 
with the 21520 cm -1 line (Fig 3.11 and 3.12). The lines at 
21543 cm - * and 21450 cm - '* are it -polarized E.D. transitions as 
they are present in then-spectrum but absent in a and axial 
spectra. They can be assigned to the two expected r— £ 
transitions. The line at 21520 cm -1 present in a and axial 
but absent in the it-polarized spectrum, is a-E.D. and can be 
attributed to the 1~1 ^3 transition. The shoulder 21518 cm” 

can not be assigned to the P[ transition as it does not 
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comply with the expected oolarization behavior. Thus ou+ 
of four, only three Stark comoonents of the D^-term are 
identified and the shoulder at 21518 also remains un- 
exolained . 

absorption 

o D 

Four fairly strong and well separated lines at 25416 cm -1 , 
25241 cm 1 , 25027 cm 1 and 24953 cm - '*' and a broad shoulder 
at 25258 cm *" associated with the 25241 cm - '*' line are present 
in the spectral region 24900 cm"'* to 25500 cm - "*" (Fig 3.9 and 
3.10). A survey of the literature and our calculations 

r p p 

indicate that the terms and must lie in the 

ranges 24300 cm -1 -24500 cm - '*', 24896-25400 cm - *' and 25900 cm - 
26100 cm - '*' respectively. Consequently, the observed spectral 
lines should correspond to the ‘ term. Theoretically a 
term with J=6 will split into 10 Stark components in the S^ 
symmetry field according to the decomoos ition (31^+4^ +3f^). 
It is reasonable to assume that — fj' transitions will not 

appear for this group because M.D. transitions are forbidden 
for Aj>+1. The line at 24953 cm - *' is very sharp and appears 
in a and axial spectra only. It can, therefore, be assigned 

to the E.D., 7 F Q ([p ^ L 6 ^ HP transition. The line at 

2502 7 cm -1 appearing only in it -polarization is it-E.D. type 

and can be assigned to the F Q ( ) L 6 (f^) transition. 

The line at 25241 cm -1 and its shoulder at 25258 cm -1 are 
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similarly o-E.D. tvoe . The relatively strono line at 
25241 cm -1 can be attributed to th^ Tp o < U > — Vrp 
transition. The shoulder is somewhat difficult to account 
for, because theoretical calculations predict a ^L^Cfp level 
just around this position, but the polarization character- 
istics of this 'shoulder' do not comply with the r i — T 
transition. The line at 25416 cm"'*' is very broad (FWHM 
34 cm ■*■) and strong, but its peak position in the it and a 
spectra differs by 3 cm - '" and the peak position in the axial 
spectrum matches with that of the o-polar i zation . In the 
n-polarization, the peak at 25418 cm - '*' can be assigned to 

the ti-E.D., 7 F q (^) ^ L 6^^2^ transition and the cr-polarized 

25415 cm - ' 1 ' peak can be attributed to the o-E.D., 7 F q ( Hp 

transition. Thus, out of the 10 expected Stark levels 
of , five are identified. These assignments are in 
qualitative agreement with the theoretically calculated 
structure of L^. However complete identifica + ion of the 

multiplet is not possible due to lack of sufficient 

experimental data and because precise calculations for the free ion 
eneraies were not possible. All observed absorption transi- 
tions are shown in Fig. 4.1. 

2. FLUORESCENCE SPECTRA 

Resonance excitation of Eu by any of the Ar laser 
lines is not possible because of the mismatch of energies. 

However LiYF 4 :Eu 3+ fluoresces when excited with most of 
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the Ar + -laser lines. The 457.9 nm (21838 cm h laser 

5 34- 

line can excite levels of Eu by phonon assisted 

processes. This excitation produces fluorescence from 

and D as well. The exact mechanism of energy transfer 
o ■ 

5 5 5 

(radiative or non-radiative ) from to D-^ and D q multi- 

plets could not be ascertained because of spectral range 
limitations of the monochromator, but radiative transitions 
are expected to be primarily responsible for this transfer. 
Similarly the 514.5 nm (19436 cm ) line can excite and 

^D 0 levels. But the fluorescence with 514.5 nm excitation 
is relatively weaker than the fluorescence with 457.9 nra 
excitation. The 5 D q level at 17271 cm -1 falls in the tun- 
ing range of the R 6G dye* Attempts to excite from the 

ground level 7 F o using dye laser were not successful because 

of severe restriction on the F q ' D q transition. The 

excitation spectrum monitored with D q F^(|^) transition 

at 16840 cm ~ 1 shows that the 5 D Q level can be excited quite 
efficiently with the dye laser lasing at 16938 cm 1 and 
17640 cm~h The excitation at 17640 cm ^ is a phonon assisted 

one while 16938 cm -1 corresponds to the '^(fp 5 D 0 (fp 

transition. We discuss below the symmetry assignments for 
various groups of fluorescence lines. 

V_ 7 f c fluorescence 

This transition is forbidden by many selection rules (ie 
(J=0 to J’=0 is forbidden for E.D. and M.D. both and Pj_ fj 
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* 

is forbidden for E.D.). This was the reason why it was 
not observed in the absorption spectrum. A very weak and 
broad line at 17271 cm corresponding to the highly for- 

bidden D (j^) ^o *1 Transition was observed when the 

sample was excited by the 457.9 nm line of the Ar + laser 
and 16938 cm ^ lines of R 6G laser. The detection of this 

line was made difficult due to the presence of a Hg line at 

—1 s 

17269 cm . Another evidence of the position of the D q 

level comes from the excitation spectrum, which shows a 
sharp excitation at 17270 cm (Section 4.3 of Chapter 
III). Polarization characteristics of this line are 
consistant with a pure o-M.D. transition. Additional weak 
but persistent lines observed at 17292 cra"\ 17250 cm - ''' and 
17231 cm - "'' can be associated with the vibronic lines (c.f. 
Table 3.8) . 

5 D 1 - 7 F. FLUORESCENCE 

I q 

This group has two lines at 19043 cm - ' 1 ' and 19021 cm - ' 1 ' 

7 5 

in agreement with the F c - absorption spectrum (Fig 3.19). 
The 19043 cm - ''' line corresponds to the ( Fji ) - 7 F Q ( ) , M.D. 
transition, while the line at 19021 cm - ' 1 ' corresponds to the 
(r^)- 7 F Q ( fp , E.D. transition. The position of both the 
Stark components of " are well established. 

5 D r - 7 F Q FLUORESCENCE 

The spectrum corresponding to this arouo of transitions 
could be obtained only with the 457.9 nm excitation. This 
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group has 4 sharp and medium strona lines at 21543 cm” 1 , 

21520 cm 1 , 21470 cm \ 21450 cm 1 and 5 verv weak lines at 
positions 21570 cm -1 , 21509 cm” 1 , 21463 cm ” 1 21458 cm -1 and 
21408 cm 1 (Fig 3.13 )• The lines 21543 cm 1 and 21450 cm ” 1 
are strong in n-spectrum and very weak in a and axial 
spectra. These are n-E.D. transitions attributed to the 
(~p(2) - an d [7,(l)-l"^ transitions. The line at 21520 cm” 1 , 
strong in o and axial spectra, but weak in the Tt-polarization, 
corresponds to the o-E.D. r 7 ^- transition. It has a small 
contribution from tc-M.D. transition as well. The 21470 cm ” 1 
line can not be assigned to the remaining Stark component 
VR ) , because it is not consistent with the other group 
of transitions. All of the weak lines are Raman or vibronic 
lines with the phonon energies 268 cm” 1 , 329 cm” 1 , 375 cm 1 , 
380 cm ” 1 and 430 cm” 1 . The line at 21470 cm ” 1 is quite 
different from the Raman lines and corresponds to the Ar + 
line at 465.8 nm. Thus all the lines of this group are 
explained. The observed fluorescence transitons terminat- 
ing on the 'F are shown in Fig. 4. 2 . The Dp(J|) 

component is not yet identified. 

5 D - V F, FLUORESCENCE 

q 

The 7 F 1 multiplet has two Stark component (1"^, f^j 4 ) • 

The observed spectrum has two sharp and extremely strong 
lines at 16938 cm" 1 , 16842 cm ” 1 and a weak shoulder at 



17271 CM 
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16933 cm 1 (Fig 3.24 ) for excitations by 457.9 mm, 514.5 m 
and R 6G (Tables 3.5, 3.6, 3.7). The line at 16938 cm” 1 is 
a ix-M.D. dominant transition with a small contribution from 
cr-E.D. and corresponds to the transition. This gives 

the position of the ^F^(Pg) level at 333 cm” 1 . The line at 
16842 cm” 1 appearing strongly in c-polari zation, but weakly in 
the axial and n-spectra, is a o-M.D. transition corresponding 
to the fi-fT transition. This would establish the second 

y — 1 

Stark component (l^) of F^ at 429 cm” . The appearance of 
the 16842 cm” 1 line in it and axial spectra is difficult to 
explain. The observed intensities are however weak. The 
shoulder at 16933 cm” 1 remains unexplained. 

5 d 1 - 7 f 1 FLUORESCENCE 

Theoretically we expect four lines for this group. We 
have observed three fairly strong lines at 18710 cm 1 , 

18687 cm” 1 and 18593 cm” 1 (Fig 3.19 and Tables 3.5; 3.6). 

The lines at 18710 cm -1 and 18593 cm” 1 are dominant cr-E.D. 
transitions with a fair contribution of tc— M.D. character. 

They correspond to the H. - q and r 3 -17 transi tions 
respectively. The line at 18687 cm” 1 is strong in ic-polar- 
ization but very we/ak in a and axial spectra, indicating 
a n-E.D. dominant character. The only possible assignment 
for this line is r 3 -g. The c -n transition has not 
been observed. These transitions re-confirm the two Stark 
components of 7 F^ at 333 cm 1 (Fg) and 429 cm 1 (!"£)• 
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i 2 “— 1 FLU0RE5CENCE 

Eight lines are observed for this group in the spectral 
range 21021 cnf 1 -21209 cm -1 (Fig 3.13 and Table 3.5). The 
line at 21209 cm’*'*' is dominantly n-M.D. and corresponds to 
the ( 2 ) — transition. The 21187 cm”'*' line is dominant 

it-E.D. but with a fair contribution of a- M.D. character. 

This line can be associated with the “[3 transition. The 
line at 21158 cm”'*' corresponds to the Ar + line at 472.7 nm. 
The fairly strong line at 21145 cm”'*' is a cr-E.D. (dominant) 
transition with a fair contribution of it-M.D. character. 
Consistent with this observation, this line could possibly 
be assigned to the ~ ^3 * P3 ~ and ” ^3 transitions. The 
last two assignments are ruled out by energy difference 

considerations. The assignment would put the “^(l"^) 

Stark level at 21478 (21145*333) cm” 1 . This assignment is 
further confirmed by the o-M.D. transition Pj_ ~ Ig appearing 
at 21048 cm ” 1 (21048+429 = 21477 cm” 1 ). Out of the remaining 
three lines, the line at 21117 cm 1 has equal intensity in 
it and a polarizations and is stronger in axial spectrum 
indicating comparable contributions from a-E.D. and it-M.D. 
transitions. This line can be associated with fT)(l) 
or with fg^) - transition. More exact assignment is not 
possible for this line. The lines at 21090 cm ” 1 and 21021 cm 
are it-M.D. and it-E.D. transitions attributed to the 
and Q(l) “ transitions respectively. All the observed 
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7 

fluorescence transitions terminating on the F-^ multiplet 
are shown in Fig 4.3. 

5 D o - 7 F 2 FLUORESCENCE 

Using dye laser excitation, we have observed three lines 

for this group at 16382 cm - ' 1 ' 16295 cm - ' 1 ' and 16122 cm"’' 1 ' (Fig 

3.25 Table 3.7). With other excitations these lines fall in the 

5 7 5 7 

overlapping regions of the Dj,- F^ and D^- F^ transitions. 

The lines at 16382 cm ^ and 16122 cm - ' 1 ' are dominantly n— 
polarized E.D. type and one assigned to the °D ( rp - 'f 2 

f^(l) and “ > D q ( ) — 7 F ? J~^(2) transitions respectively. This 

m 7 

assianment would place the two Stark levels of F^ at 
889 cm - ' 1 ' and 1151 cm - ' 1 '. The line at 16295 cm -1 is a dominant 
tc-M.D. with comparable contribution of cr-E.D. character and 
would be consistent with the f"^ - assignment. We infer 
the position of the component of F? at 976 cm - . The 
- j"^ transition could not be seen. 

5 D 1 - 7 F 2 FLUORESCENCE 

The seven lines observed in the spectral range 17840 cm 
18155 cm - ' 1 ' belong to this group (Fig 3.20 ; Tables 3.5 and 
3.6). The weak line at 18142 cm -1 however does not belong 
to Eu^ + . The very weak but pure it-E.D. transitions at 
18153 cm -1 and 17892 cm -1 are the -f^(l) and ri - ^2^ 2 ^ 
transitions terminating at 890 cm ^ and 1152 cm ^ respectively® 
The lines at 18132 cm -1 , 18067 cm -1 and 17869 cm -1 are dominant 
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tc— M.D. with some contribution of g-E.D. character and are 

assigned to the f" 7 ^ -r 2 ( 1 ), fi -r 3 and r 3 -r 2 ( 2 ) transitions . 

The line at 18044 cm is dominantly c— M.D. and is attributed 
to the r 3 -I 3 transition. The line at 17846 cm " 1 is the 
ti-m.d. (r 3 -n ) transition and gives the position of the 
7 Fi(q) level as 1175 cm" 1 . 

" >D 2~ ?F 2 FLUORESCENCE 

Theoretically, we expect 16 lines for this group. We 
have observed 13 lines in the spectral range 20270 cm " 1 - 
20650 cm 1 (Fig 3. 14, Table 3.5). The line at 20629 cm " 1 
with mixed o-E .D. (dominant ) and -rr-M.D. character can be 
assigned to the -f^(l) transition and the 20589 cm " 1 line 
with 'rc-E.D. character is assigned to the -fj(l) transi- 
tion. The lines at 20567 cm" 1 , 20502. cm " 1 and ,20474 cm " 1 
are cr-E.D. (dominant) with some tc-M.D. character, while 
the line at 20545 cm " 1 is dominant -Jt-E.D. with a small 
contribution from o-M.D. They respectively correspond to 

the r 2 <2) -rj> ir-i3.fa (i) -fi and r 3 -13 transitions and 
give the position of the terminal level ^pd"^) at 976 cm" 1 , 
977 cm" 1 , 976 cm" 1 and 975 cm" 1 respectively. The 20390 cm" 1 

line is extremely weak and can be assigned to the 1^(2) - 

7 

f^(2) cr-M.D. transition suggesting the terminal level 
g< 2 > to lie at 1153 cm" 1 . The broad line at 20369 cm"" 1 
unpolarized and may correspond to the two 


appears 
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transition f 3 -[*(2) and ( 2 ) - f] terminating on the 1151 cm" 1 

and 1174 cm" 1 levels of 7 F 2 respectively. The weak lines 
at 20344 cm 1 and 20271 are assigned to the (a E.D.,n: M.D.) 
r 3 -fj and it-E.D. r 2 (l)-fjf transitions respectively giving 
the positions of the terminal levels at 1176 and 1179 cm" 1 al- 
though the polar ization behavior of the 20271 cm" 1 line is not 
quite consistent with this assignment. Out of the remaining 
three lines, the one at 20493 cm" 1 being quite sharp and stroncr 
corresponds to the 488.0 nm of Ar + line and the extremely 
weak lines at 20578 cm" 1 and 20455 cm" 1 remain unexplained. 

Pig. 4.4 shows all the observed transitions terminating 

7 

on the F 2 multiplet. 

5 d ,- 7 f 3 FLUORESCENCE 

The J=3 multiplet splits into five Stark components in 
S^ symmetry (^+2^ +2 4 ) . We have observed 4 well sep- 

arated lines for this group (Fig 3.26 , Tables 3.5, 3.6 and 
3.7) in the spectral range 15300 cm" 1 - 15420 cm" 1 . The 
line at 15411 cm" 1 is tc-E.D. dominant and corresponds to 
the [^-{^(l) transition. It assigns the level at 1860 cm" 1 
to one of the components of 7 F 2 « The T5398 cm ^ an ^ 

15317 cm" 1 lines with dominant o-E.D. character with fair 
contribution from tc-M.D. is assigned to the transi- 

tions. From these assignments, we deduce the f"* 3 (l) and 
r 3 (2) components of 7 F 2 at 1873 cm" 1 and 1954 cm" 1 respectively 



479 



(♦Indicates wavenumbers corresponding to two 
transitions ) 
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The line at 15369 cm ^ is a pure a-M.D. (iy ~ pp transition 
and gives the position of the 'F^(P^) level as 1902 cm . 

V 7f 3 FLUORESCENCE 

Out of the expected ten lines belonging to this group, 
we have observed nine in the spectral range 16950 cm*"'*' to 
17200 cm" 1 . There are 4 lines of mixed nature (o-E.D. + 
it-M.D.) with cr-E.D. dominant at 17172 cm” 1 , 17162 cm \ 

17119 cm -1 and 17092 cm -1 . These lines are attributed to 

the r i r 3 -fjCl), r 3 -Fi and V[ - 1 ^( 2 ) transitions 

7 

respectively qiving the terminal levels of the F^ multiplet 

at 1871 cm -1 , 1859 cm -1 , 1902 cm -1 and 1951 cm -1 . The line 

at 17185 cm -1 , overall weak but dominantly it -E.D. corresponds 

to the -p, (l) transition terminating at 1858 cm h The 

three lines at 17149 cm -1 , 17069 cm -1 and 17006 cm " 1 are 

predominantly te— E.D. with small intensity in a— polariza— 

tion correspond to the P 3 -[" 3 ( 1 ), F 3 “ 1 ^( 2 ) and 

transitions respectively. The terminal levels are 1872 cm , 

1952 cm " 1 and 2038 cm" 1 . The remaining one line has 

comparable contribution from 71 — M.D. and 0 — E.D. transition 

and is attributed to the f 3 -1^(2) transition with the 

7p ( p ) at 2038 cm -1 . Thus this group gives all the Stark 
3 2 

components of 7 F 3 and confirms the assignments of the 

5 V 7f 3 gr0Up ' 
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assigned to the P 3 -r 3 ( 2 ) transition, terminating at 
1954 cm "*■ . The 19548 cm ^ is dominantly n-E.D. f^Cl) - transi 
tion, terminating at the 1902 cm -1 level, but its weak 
appearance in the a and axial spectra is not consistent with 
this assignment. The lines at 19525 cm ” 1 and 19480 cm ” 1 
are dominantly tc-M.D. assigned to the P^ -f^(2) and p^ -pCo), 
transitions, terminating at 1954 cm -1 and 2040 cm -1 respect- 
ively. The transition corresponding to 19496 cm -1 is n-M.D. 
with comparable contribution from cr-E.D. and is attributed 
to the r 2 (l) -[ 3 ( 2 ), transition with P^ 1954 cm - " 1 . The 
last line at 19439 cm -1 is extremely weak but pure n-E.D. 

<ri ) transition, terminating at 2040 cm”' 1 . These 

assignments are in agreement with the assignments made in 

the previous two groups. Fig. 4.5 shows all transitions 

7 

terminating on the F^ multiplet. 

5 D c - 7 F 4 FLUORESCENCE 

The multiplet with J=4 splits into 7 Stark components 
in a crystal field of the S^ symmetry (31^+21^+21^). We 
have observed 4 stronq peaks and two very weak lines for 
this group in the soectral region 14200 cm ^-14700 cm 1 
(Fig 3.27 and Tables 3.5, 3.6 and 3.7). The spectra in 
the range 14500 cm -1 - 14700 cm ” 1 are not included in the 
Fig 3. 27 • The two weak lines at 14664 cm -1 and 14400 cm 1 
are cr-M.D. - 1"^ transitions, terminating on the (^(l) and 
p[(2) levels at 2607 cm -1 and 2871 cm -1 respectively. The 



G<2) 
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lines at 14458 cm ^ and 14259 cm -1 are (tc-M.D.+ a-E.D.) 
n transitions. They predict the two p 3 components of 

7 _i 

at 2813 cm and 3012 cm . The remaining two lines 
at 14366 cm ^ and 14294 cm""'*' are tc-E.D., I~1 -II transitions 
and assigned the positions of the two components at 
2905 cm - "*" and 2977 cm - '*’ . 

5 d 1 - 7 f 4 FLUORESCENCE 

We have observed 11 lines belonaing to this group out 

of 14 expected lines (Fia 3.22 Table 3. 5 and 3. 6 ). The 

line 16436 cm - '*' observed only with the 457.9 cm 1 excitation 

is extremely weak a-M.D. type attributed to the r, -qu) 

transition yielding 7 F 4 (f^) level at 2607 cm 1 . The line at 

16413 cm"" 1 is dominantly o-E.D., r^-fj(l) transition with 

terminal level n at 2608 cm~^-. The line at 16229 cm 1 

having comparable contribution from rc-M.D. and o-E.D. 

corresponds to the Tj, -f^(l) transition, terminating at 

2814 cm -1 . The line at 16208 cm -1 , it-E.D. (dominant) + er-M.D. 

7 t mm f 

is the r*3 -g (l) transition suggesting the F 4 (r 3 ) level 
at 2813 cm -1 . The lines at 16149 cnT 1 and 16138 cm" 1 are not 
well resolved, but their observed polarization behavior 
predicts them to be (o-E . D.+it-M .D.) -1^(2), and it-E.D. 

(R l)) transition corresponding to the F 4 (P[) at 

2872 cm" 1 and 7 F 4 ([p at 2905 cm" 1 respectively. The 
extremely weak line at 16064 cm -1 shows up only in the it- 
polarization and hence can be assigned to the -f^(2), 
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it-E.D. transition. The strong line at 16042 cm 1 has 

comparable contribution from er-E.D. and it-M.D., r 3 -rj(a) 

transitions. Both the transitions assign to the {“2,(2) level a 

position of 2979 cm The line 16009 cm 1 is very s trono 

and dominantly tc-E.D. type assigned to the j"^ - JJ(?) transi- 

tion with 3012 cm as the position of F 4 ( f^) line. A 

shoulder at 16032 cm -1 associated with the 16042 cm -1 line 

also corresponds to the 2 ) transition and gives the 

r,(2) level a position of 3011 cm -1 . The line at 15940 cm *" 1 

is (ix-M.D.+cr-E .D. ) type and can be assioned to the {"3 -f^(3), 

transition givina rr<3) a position of 3081 cm"* 1 * . Thus all 

7 

the seven Stark comoonents of the are identified. 

5 D^- 7 F / , FLUORESCENCE 

q. 

The spectrum in the spectral range corresponding to 

this group is somewhat complex. The strong transitions 

5 7 

belonging to the D, - F^ may mask many of hhe lines of the 

c n 

D n - F. grouo. Some of the lines are not well resolved, 

2 4 

others are very weak. We have been able to identify some 

16 lines belonging to this group out of the 28 lines 

expected theoretically. There are many extra lines which 

3+ 

could not be explained in terms of Eu transitions and 
remain unassigned (Fig 3.16 , Table 3.5). There are two 
extremely weak lines at 18912 cm 1 and 18842 cm 1 (not 
shown in the Fig), which correspond to the -1^(1) ancJ 
P 2 (l) -f^(l) transitions if the energy differences are 
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matched. The terminal level p[(l) is at 3608 cm -1 . The 

two weak lines at 18666 cm -1 , 18636 cm -1 are (tc-M.D. 

( dominant ) + a-E . D . ) and a-E.D. respectively and are assianed 

to the o -m 1 ) and r 2 u) -qc l) transitions re socc tive ly . 

The line at 18619 cm 1 appearing fairly strona in a-polar- 

ization is somewhat masked by the 18595 cm -1 line of the 
5 7 

Df- F^ group, corresoonds to the 9s -G< l) transition with 
7 F 4 (r 2 ) level at 2901 cm _1 . The line at 18574 cm - ^ a clear 
shoulder in the m-and a-polari zation but masked in the axial 
spectra by the 18595 cm 1 strong line may correspond to the 
n -f2 (i) transition with[7(l) at 2905 cm Again at 18565 cm ' 
isaweak shoulder appearing more clearly in a-polari zation, which 

y 

corresponds to the r^(2) transition with F^ 

level at 2978 cm -1 . The line at 18542 cm -1 is a-E.D. (domi- 
nant J+m-M . D . 9 -F< 2) transition with r 2 ( 2) level at 2978 
cm -1 . The line at 18528 cm - "'' is tt-M.D . (dominant)+a-E .D. , 

j 

r 2 (2) -q (2) transition and confirms the position 

at 3015 cm . By energy difference consideration alone, 
the line at 18505 cm -1 can be assigned to the -f^(2) 
transition though polarization behavior does not quite 
satisfy this assignment. There are two overlapping lines 
at 18471 cm -1 and 18467 cm” 1 which correspond to the ("^(l)- 
r 2 a> and 0 -8 (2) transitions. The line at 18436 cm 1 
■rc-M.D. (dominant )+a-E . D. is a superposition of the r 3 -90) 
and r,u) -qc 2) transitions. The line at 18394 cm 1 
corresponds to the q -9(3) transition according to energy 
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difference consideration alone . All the a s s icf nine nt s in 

this group support the assignments made for the 5 D - 7 F„ and 

o 4 

5 7 

^4 groups. Fia 4*6 shows observed fluorescence transi- 

7 

tion terminatina on the ? 4 multiple t* 

5 d - 7 f c fluorescence 

o 5 

7 

The F^ multiplet splits into 8 Stark components in S^ 
symme try ( 3r 1 + 2g + 3m 4 ) • We have observed only four lines 
in the spectral ranae 13200 cm -1 - 13500 cm -1 (Fig 3.28 , 
Table 3.7). These lines cou.ld bo observed only with R- 6 G 
dye laser excitation at 16938 cm - '*' . The line at 13485 cm - " 1 " 
is dominantly it-E.D. type and is attributed to the F^ -F,(l) 
transition. The two lines at 13472 cm - '*" and 13271 cm ^ are 
o-E.D. types and are assigned to the -1^(1) and F^ -1^(2) 
transitions respectively. The line at 13^19 cm ^ seen in 
cr-polari zation only but extremelv weak can be assigned to 
the Pj, -(" 3 ( 3 ) transition. These four transitions give the 
4 Stark levels, F>(l), F^Cl), 1^(2) and f^(3) of the 7 F 5 
multiplet at 3786 cm ^ , 3799 cm ^ , 4000 cm and 4052 cm 
respectively . 

5 d 1 - 7 f 5 FLU ORESCENCE 

We have observed 15 lines in the spectral range 14900 
cm -1 - 15300 cm -1 (Fig 3.23 and Table 3.5). This group 
could be recorded only with 457.9 nm excitation. The line 
at 15283 cm -1 could not be identified. The line at 15245 cm 
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is cr-E.D. (dominant )+ tc-M. D . type and is attributed to the 
ri -ij< l) transition with the ["^(l) level predicted at 
3798 cm'*' 1 '. The 15236 cm”' 1 ' line appears to be unpolarized 
and energy wise can correspond to the -f^( 2 ) transition 
belonging to 5 D o - 7 F 3 group or the P 3 -^(l) transition 

c n _ 1 

belonging to ^D^- group. The line at 15223 cm is 
'rc-E.D. type and is assigned to the P 3 -f^(l) transition. 

The lines at 15214, 15203 and 15054 cm ” 1 are extremely weak 
and could not be assigned to Eu 3+ transitions. The line at 
15191 cm -1 predicts a level of 7 F^ at 3852 cm 1 or a 

level at 3830 cm” 1 , both are not supported by any other 
evidence and conseguently this line is left unassigned. The 
15177 cm ” 1 line is cr-E.D. and corresponds to the P 3 -fjf(l) 
transition. This assigns the F^(l) level of the Fy Mul- 
tiple t at 3844 cm” 1 . The extremely we ak , cr-E . D .-type line at 
15044 cm ” 1 is assigned to the [J- fj( 2 ) transition, and gives thef^( 2 ) 
line of 7 F 5 , a position of 4000 cm” 1 . The relatively 
strong line at 15022 cm ” 1 (cr-E .D ,+cx-M.D . ) is attributed to 
the r 3 -{ 3 ( 2 ) transition. The lines at 15014 cm ” 1 and 
14994 cm ” 1 are cr-E.D. transitions assigned to the Pj -1^(2) 

and P, -P(3) transitions respectively. The assignments 
X O 

predict the positions of the f^( 2 ) and P^ 3 ) levels of F^ 
at 4007 cm ” 1 and 4050 cm ” 1 respectively. The 14971 cm 
line is dominant tc-E.D. type with some contribution from 
a-M.D. and corresponds to the P 3 -T^ 3 ) transition. 
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5d 2~ ?F 5 flu ORESCENCE 

This group of fluorescence lines is observed in the 
spectral range of 17370 cm” 1 - 17750 cm -1 (Fig 3.17 , Table 
3.5). The Following lines are identified as belonging 
to this group. The lines at 17745 and 17719 cm” 1 are 
(it-M.D.+o-E .D. ) and (c-M. D.+it-E .D. ) types corresponding 
to the f 2 (2) -rj(D and r 3 -r 3 (l) transitions respectively. 

The 17691 cm 7 line is a tc-E.D. type and can correspond to 
the r x -f^d) and p 2 (2) -rj( l) transitions. It is not 
possible to make a more definite assianment for this line. 

The line at 17678 cm -1 is (n-M. D.+cr-E . D . ) and is assigned 
to the P^ -rj(l) transition. The three lines at 17605 cm - "!; 
17533 cm -F and 17443 cm - ' 1 " are it-E.D. type and are assigned 
to the P 2 (l) -f^(l), r 2 (2) -fj(2), fj^(l) “fid) transitions 
respectively. Further from cr-E.D. transitions lines at 17480 cm 
17451 cm - '*' , 17428 cm - ' 1 ' and 17400 cm - ' 1 ', are attributed to 

the Px -r 3 ( 2), r 2 ( 1) -fj(2), P x -r 3 (3) and fyl) -JJ( 3 ) 

transitions. Thus out of the eight Stark components of 

7 F C we are able to qive tentative assignments to six. The 
5 

7 

observed fluorescence transitions to the F^ multiplet are 
shown in Fig 4.7. 

5 D 2 ~ 7 F C F LUORESCENCE 

5 5 

We could not observe any fluorescence from D Q and to 
the 7 F^ mu It ip let. The 5 D 2 ~ 7 F 6 fluorescence has been 
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observed (Fig 3.18), but could not be analyzed. These lines 
are closely spaced, broad and not well resolved. 

Table 4.1 lists the mean positions of the energy levels, 
the tentative symmetry assignments and the CG's of the 
multiple ts as obtained from experimental results. The 
procedure employed for obtaining these results is outlined 
in appendix III. 

EXTRA-LINES 

The fluorescence spectra recorded from our sample have 

3+ 

many lines which apparently do not belong to Eu . These 
lines are marked a, b or c in Figures 3.13 to 3.28 and 
Tables 3.5 to 3.7. The lines marked ’b' are Hg-lines due 
to the room light while those marked ' c' correspond to Ar + 
lines. There is a very strong line at 20866 cm - "*' and a weak 
line at 20788 cm - , which may correspond to the P Q (rp - 
3 H 4 (rp and ^ P 0 (rp ~ ^ 4 (^ 3 ) transitions of Pr 3+ . This 
indicates the presence of small Pr 3+ impurity in the crystal 
under study. Other extra lines which may originate from 
Pr 3+ impurity are 20706 cm - ^ (vibronic), 20671 cm ^ vibronic 
and 18580 cm -1 ( 3 P o ~ 3 H 5 ). The remaining extra lines marked 
'a' have not been accounted for in this work. 

3. CALCULATION OF ENERGY LEVELS 

As mentioned earlier, H c ry <H s- 0 +H Coul but H cry is not 
so small to be confidently treated within the framework of 
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TABLE 4.1 

OBSERVED POSITIONS OF ENERGY LEVELS OF LlYF^Eu 3 * 


S.No. 


Multiplet 


Symme try 
Assignment 
( tentative ) 


Energy C.G. 

(cm - *) (cm*"*) 



T 

T 

r r 

Ig 

I? 

r? 

GU 

IT 

f?,4 

I? 

IT 

17,4 

I? 

I? 

17 


0 

429 

333 

1175 

1152 

976 

889 

2038 

1954 

1902 

1873 

1859 

3081 

3012 

2979 

2905 

2871 


0 

365 

1033.6 


1921.9 


2899.3 
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n, 4 

2813 




IT 

2608 


6. 

\ 

rr 

- 

3948 



t?.4 

4052 




rr 

- 




rr 

4008 




(T,4 

4000 




rr 

3844 




ta.4 

3799 




rr 

3786 


7. 

\ 

rr 

17271 

17271 

8. 

\ 

r? 

19043 

19028.3 



t?,4 

19021 


9. 


rr 

21543 

21502.4 



f^4 

21520 




rr 

21479 




rr 

21450 


10. 

% 

rr 

25418 

* 



ir.4 

25415 




rr,+ 

25241 




rr 

25027 




rr.+ 

24953 


♦Two 

and three 

fl levels of 5 l^ 

could not be 

obtained . 
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the first order perturbation. However the exact calculation 
requiring simultaneous treatment of these interactions is 
not only cumbersome but has its own shortcomings; for example 
the simultaneous parametrization of the crystal field and 
free ion interactions tends to accommodate the inadequacies of 
each other. Although such an approach has been used in some 
cases but the usual procedure is to carry out such calcula- 
tions in two steps. In first step the spin-orbit and the 
residual Coulomb interactions are considered and the wave 
functions so generated are used to obtain the crystal field 
splittings . 

3.1 FREE ION CALCULATIO NS 

The following procedure was adopted to calculate the 
positions of the free ion energy levels. The free ion energy 
matrices for the Eu^ + (4f^ configuration) were constructed 
in the intermediate couplinq scheme following the procedure 
outlined in appendix III. The free ion Hamiltonian included 
the Coulomb, the spin-orbit and the configuration interactions. 
The angular parts of the Coulomb interaction were directly 
taken from Ref. [8] of Chapter II, while -angular parts of the 
spin-orbit and the configuration interactions were calculated 
using the expression^ discussed in Chapter II. We have no 
experimental data for multiplets with J more than six; conse- 
quently the programs were written for multiplets with J = 0 
to 6 only. These programs use the subroutine 'EIGEN* obtained 
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from 'IBM 360 Scientific Subroutine Package' to diagonalize 
the energy matrices and to arrange the eigen values in the 
ascending order. The free ion matrices were constructed 
one for each J value. Since the sizes of the free ion energy 
matrices were large requiring a huge amount of data to be 
fed to the computer to build up these matrices, it was desir- 
able to check the correctness of the programming process. 

This was done by comparing some of our calculations with 
those applicable to Eu^ + in aqueous solutions [l]. There was 
almost complete agreement in general, however, a few discrep- 

R R CL 

ancies exist. Our calculated positions of D 4 , F 4 and I 4 
levels are 27672, 33672 and 33893 cm~^ whereas they are placed 
at 27670, 33651 and 33914 cm ^ respectively in Ref.[l]. We 
checked our programs and energy matrices several times, but 
could not find any reason for this discrepancy. 

Experimentally, the free ion energy levels are obtained 
from the centres of gravity (CG's) of the observed positions 
of the Stark components of a given multiple t. Usually some 
components are missing in the spectra either because the 
relevant transition may be forbidden or only weakly allowed. 

In such cases the approximate positions of the missing levels 
can be estimated on the basis of their calculated positions. 
The approximate CG's so obtained can be further improved with 
each refinement in the crystal field calculations. The 
results presented in this section are obtained by considering 
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only those multiplets (except 7 F 5 ) for which all the Stark 
components could be identified. The free ion level posi- 
tions of 7 F o , 7 F 1 , 7 F 2 , 7 F 3 , 7 F 4 , 5 D q , 5 D 1 and 5 D 2 are well 

established from our experimental results. The missing two 
1 

levels of ^(["^(3) anc | 1 ^( 2 )) were taken at their calculated 
positions in reference [ 2 ] . A program for free ion para- 
meterization can be written using the minimization subroutine 
'MINFUN 1 (described in appendix III). The large size of the 
free ion matrices, however, makes any minimization exercise 
impractical (and expensive) for Eu 3+ , as it requires enor- 
mously large CPU time. We adopted the following approximate 
procedure to reduce the discrepancy between the observed and 
calculated energies of free ion levels. We started with the 
free ion parameters given in reference [ 2 ]. The starting 
parameters gave an r.m.s. deviation of nearly 55 cm The 

free ion parameters were changed one at a time in small 
steps and the r.m.s. deviation calculated at each step. Thus 
for each parameter, the minimum r.m.s. deviation and the 
corresponding values of the parameters were obtained, while 
rest of the parameters were kept at their starting values. 
When these new parameters were used for free ion energy 
calculations, the r.m.s. deviation first went-up to 96 cm 
However, slightly modified procedure of varying the para- 
meters brought the r.m.s. deviation down to 24 cm -1 . Some- 
what different procedures gave slightly different r.m.s. 
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deviation. All these trials indicate that we do not have 
a unique set of parameters but the parameters obtained in 
this manner (Table 4.2) are not too far off the mark. These 
parameters were used to obtain the zero-order wavef unctions 
needed to carry out the crystal field calculations. Table 
4.3 gives the free ion parameters, free ion energies and 
eigenfunctions . 

3.2 CRYSTAL FIELD CALCULATIONS 

(A) CALCULATIONS WITHOUT J-J MIXING 

In order to obtain crystal field split energy levels, 
the crystal field energy matrices were constructed for the 
three representations ir rj. 17,4 for each multiplet separ- 
ately. This procedure neglects J-J mixing effects. The 
zero-order wavef unctions obtained in the previous section 
were used for this purpose. The angular parts of the crystal 
field matrix elements were calculated using the expressions 
given in Chapter II. A general program was written to 
calculate these angular parts by defining two quantities 

ATOM(k) = (-l) S+L,+J+k V[(2J+1)(2J'+1)] ^ g j* 

x <f 6 uwslI |U (k) | If 6 U'W' SL ’> 1 

)<f | |C^ I |f> 

J' 


and 


F(k,q) = (-D J - J z ( 


- J z q 


2 
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TABLE 4.2 

EXPERIMENTAL AND CALCULATED FREE ION LEVEL-POSITIONS 


S.No . 

Multiple t 

Observed 
pos ition 

(cm -1 ) 

Calculated 8 
pos ition 

(cm”*) 

Calculated 13 

position 

(cm”*) 

1. 

7 f 

0 

0 

0 

0 

2. 

?F 1 

365 

369 

379 

3. 

?p 2 

1034 

1071 

1043 

4. 

?P 3 

1922 

1854 

1899 

5. 

7p 4 

2899 

2816 

2876 

6. 

?F 5 

3948 

3846 

3922 

7. 

5 °o 

17271 

17234 

17229 

8. 

\ 

19028 

19001 

19041 

9. 


21502 

21446 

21535 



r.m.s . 
deviation 

56 cm”* 

24 cm”* 


Calculated using parameters of reference [2]. 

E* = 5544 cm"”* E^ = 24.83 cm”* = 585 cm * 

= 1285 cm”* a = 20 cm”* p = -640 cm”* y = 1750 cm 


Calculated using our values. 

E * - 5549.2 cm”* E^ = 24.8 cm”* E^ = 585.2 cm 
C, = 1307 cm”* a = 17.1 cm”* p = -639.4 cm”* 

y = 1749.8 cm”* 


r 


3 


2 

i 
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The quantities <| |U^^| |> were obtained from the tables 
of Nielson and Koster, while 3-J and 6— J symbols were cal- 
culated using appropriate subroutines. Some matrix elements 
were calculated manually in order to check the computer 
programs written for this purpose. The energy matrices 
obtained in this manner were diagonalized using the ’EIGEN' 
subroutine. The optimization subroutine 'MINFUN' was used 
to obtain the best set of parameters. A brief note on the 
•MINFUN' subroutine and programming is given in appendix III. 


Out of the seven parameters (B^ ’ s ) appropriate to 
S 4 symmetry group, the axial parameters (q=0) are indepen- 
dent of the choice of x and y axes of coordinates. The 
values of the remaining four parameters (with q 4 0) are 
dependent on the choice of x and y axes. Although the choice 
of the z-axis is dictated by the crystal symmetry axis (c- 
axis), no unique choice exists for the x and y axes. A 
rotation of the coordinate system through an angle (j) about 
the z-axis leaves the axial parameters unchanged, but the 
non axial parameters transform as [3], 


[Bq] = Cos(q,<t>) [B*j + Sin (q (|>) 


Where [B k ] are the parameters in the new coordinate system 
and are those in the old coordinate system. A proper 

V_ ^ J r k T 

choice of angle (p can make one of the LB^J parameters zero, 
so that we need to determine only six parameters. In our 
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4 / '4s 

calculations presented here, we have set B_ 4 (= ) to zero* 

At one stage we tried to consider all the seven parameters 

tg. 

and also the six parameters with B 4 = 0, but the results 
were unchanged. 

The experimental values of the crystal field splittings 
were obtained by taking differences of the Stark components 
of a given multiplet with respect to its lowest Stark com- 
ponent. In this manner, a total of 23 such energy differences 
were obtained, 19 belongino to the Fj(J=0 to 5) and 4 to 
the ^Dj(j=0 to 2) multiplets. 

The multiplets with J=0 are non-degenerate and do not 
play any role in the crystal field parameterization as far 
as calculations without J-J mixing are concerned. We did 

5 

not consider the multiplet for lack of sufficient experi- 

mental data. 

Gorller et al [2], have reported results of crystal field 
parameterization for the Stark levels of the Fj(J=0 to 5) 
multiplets and obtained an r.m.s. deviation of about 3 cm”' 1 ' 
considering J-J mixing within the 7 Fj multiplets. However 
when we used their parameters to calculate the crystal field 
splittings without J-J mixing effects, the resulting r.m.s. 
deviation was approximately 23 cm and when and 

multiplets were also included the r.m.s. deviation went up 
to about 25 cm“*. We carried out a least squares fitting to 
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» A 

the six crystal field parameters (B^ was arbitrarily set 
to zero) in an attempt to bring better agreement in our 
calculations and experimental results. The subroutine 
’MINFUN’ was initialized for the ’search mode' with para- 
meters of ref [ 2 ]. This gave a large number of approxi- 
mate minima and the corresponding parameters. When these 
sets of parameters, were used as starting parameters in the 
’convergence mode’ of the 'MINFUN * , they all converged to 
values given in Table 4.4. The r.m.s. deviations^ obtained 
from our least square fittings to 19 energy differences 
(corresponding to Fj, J=0 to 5) and 23 energy differences 

e "i 

(including 3 Dj, J=0 to 2) were 17.5 and 20.1 cm"” respect- 
ively. The calculated energy splittings along with the 
experimental values and the crystal field wave-vectors are 
$ The r.m.s. deviation has been defined as 

a = {ii Aj)/n 
i 1 

Where A ^ is the differences between calculated and experi- 
mentally observed splittings; n is the number of energy 
differences used in the least squares fit. Here a point is 
worth noting that if we take n as the number of levels 
involved the r.m.s. deviation will get reduced because the 
lowest components of multiplets taken as reference to 
calculate splittings will not contribute anything to the 
numerator while denominator will have a higher value. 
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TABLE 4.4 

CRYSTAL FIELD PARAMETERS (in cm" 1 ) -WITHOUT J-J MIXING 



Set A 

Fit with 19 energy 
differences 

Set B 

Fit with 23 energy 
differences 

B 2 

o 

380.60+ (2.29) 

370.54+ (2 .23) 

< 

-818.25+ (3 .79) 

-821. 25+ (3. 78) 

< 

- 81.618^(3.79) 

- 81.452+ (3 . 79) 

B 4 

-972. 95+ (2. 98) 

-1001.0+ (2 . 92) 

B 4 

-825.72+ (2.66) 

-816 .02+ (2. 56) 

CQ 

0.0 

0.0 

b '6 

B 4 

300.73+ (7 . 54) 

272.46+ (8.41) 

r.m.s. 

deviation 

17.5 cm" 1 

20.1 cm” 1 


. In set A only 7 Fj multiple ts are considered whereas 

5 

parameters of set B are obtained when Dj multiplets are 
also included. 

. Calculated splittings using parameters of set A and set B 
differ to an extent of +2 or 3 cm -1 for some of the levels. 
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given in Table 4.5« To check reliability of our minimiza- 
tion procedure we repeated the above calculations with a 
totally random set of starting parameters. The final values 

2 4 6 

of B o , B q and B q remained the same, but some of the non- 
axial parameters changed sign. The change in the sign of 
the non-axial parameters (b£, q^O) did not modify the r.m.s. 
deviation at all. The theoretically calculated order of 
energy levels is in agreement with the experimental results. 

The relatively large r.m.s. deviation here is due to the 
neglect of J-J mixing effects and partly due to the fact 
that we have more data to account for. To estimate contribu- 
tion from J-J mixing effects we have attempted a limited J-J 
mixing calculation. 

B. CALCULATIONS WITH J-J MIXING • 

A detailed consideration of J-J mixing effects in 
3+ 

LiYF 4 :Eu needs to solve the energy matrices of order upto 
761x761, even if the matrices are broken into the three sub- 
matrices for each representation. This is a rather difficult 
job. We have carried out a limited J-J mixinq calculation 
involving only two lowest multiplets for each J from J=0 to 
J=6. Thus six crystal field energy matrices corresponding 
to three irreducible representations and two multiplets were 
constructed considering matrix elements between |JJ 2 > and 
jj’, J' z > bases with J' ^ J as well. The procedure for 
calculating matrix elements is outlined briefly in appendix III. 
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The splittings are given «itn respect to tne lowest Stark: level (observed) of a. given multioiet 

The 516 multioiet is not included in tge crVstai field least squares analysis 

The numbers in the oarenthesis" ( ) " represent respectively the real and imaginary Parts of a 
coefficient. 
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The expressions and general procedure outlined where J-J 
mixing effects are neglected are applicable here as well. 

The resulting sizes of matrices did not exceed 13x13, but 
the size of intermediate computations increased immensely. 
Therefore the programs were written to calculate and use the 
ATOM(k) and F(k,q) in the program itself and to obtain 
angular parts of the matrix elements only in the final stage, 
decreasing storage requirements. The matrices were diag- 
onalized in a similar fashion as described in the previous 
section. To take account of free ion splittings the 
calculated free ion eigen values were added to the diagonal 
matrix elements appropriately. For least squares fitting 
the comparison between observed and calculated crystal field 
splittings was made taking lowest Stark component of each 
multiple t as reference. 

When the parameters of ref [ 2 ] were used for these cal- 
culations, resulting r.m.s. deviation was approximately 

*7 

5 cm - '*' for 19 energy differences corresponding to Fj(J=0 
to 5) and when 5 Dj(J= 0 to 2) were also included it went up 
to about 14 cm -1 for 23 energy differences. The six para- 
meter least squares fits were attempted for 19 energy 
differences and then for 23 energy differences. We could 
achieve an r.m.s. deviation as low as about 2.9 cm 1 and 
13.1 cm -1 for 19 and 23 energy differences respectively. 
Table 4.6 gives the crystal field parameters obtained from 
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T ABLE 4.6 

CRYSTAL FIELD PARAMETERS (in cm" 1 ) (WITH J-.J MIXING) 



Set C 

Fit with 19 
energy differences 

Set D 

Fit with 23 
energy differences 

B 2 

367 .40+ (2.31 ) 

360.51+ ( .23) 

0 



B o 

-721.79+ (3.71) 

-723. 18+ ( .49) 



B 6 

o 

-43.919+ (4.60) 

-42. 617+ (.43) 

B t 

-927.66+ (2.74) 

-929. 6 6+ ( .12) 


-809.77+ (2.59) 

-808.83+ ( .13) 

*4 

B 4 

0.0 

0.0 

b ; 6 

237.84+( 9.00) 

238 .39+ ( .45) 

r .m.s . 
deviation 

2.9 cm- 1 

13.1 cm”* 1 


In set C only 7 Fj multiplets are considered whereas para- 

5 

meters of set D are obtained when Dj multiplets are also 
included. 

Calculated splittings using parameters of 
differ at most by + 2 or 3 cm 1 . 


sets C and D 
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these fits and Table 4.7 gives the experimental splittings, 

the calculated splittings and the crystal field wavef unctions 

(with J-J mixing effects)* We notice considerably improved 

agreement for the low lying levels (i.e. 7 Fj multiplets), 

thus establishing a significant contribution of J-J mixing 

effects for these closely spaced multiplets. The neglect 

of higher multiplets in these calculations is justified by 

relatively pure character of F wavef unctions (96'/ to 99/), 

Other workers [4,5] have also used somewhat similar approach 

3+ 

to calculate energy levels of Eu in other hosts. The 
agreement for higher levels ( Dj multiplets) is not so good 
even when J-J mixing effects are considered. A similar 
situation is encountered elsewhere [5] where somewhat 

different truncation procedure was used. The mixing of 

3 

the P states was suggested as a possible cause for this 
disagreement. This however may not be so as we have used 
intermediate wavefunctions which take care of these effects 
to a certain extent. Though we do not have any satisfactory 
explanation for this poor agreement, we nevertheless feel 
that a more exact J-J mixing calculations involving other 
nearby multiplets should improve the situation but computa- 
tional complexities have prevented us from performing such 
calculations. It should be pointed out that we accounted 
for the free ion contributions in the J-J mixing calcula- 
tions by adding to the diagonal matrix elements, the calculated 
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free ion positions but in principle one should diagonalize 
the free ion and crystal field energy matrices simultaneously. 
A comparative look of Tables 4.5 and 4.7 shows clearly that 

7 

the effect of J-mixing is prominent on the Fj multiplets 
s 

while the Dj multiplets are somewhat unaffected by this. 

This is rather strange. 

In conclusion we can say that our limited J-J mixing 
calculations using intermediate wavefunc tions seem to 

7 

adequately account for the observed splittings of the Fj 

5 

multiplets but not guite so for the Dj and other multiplets. 
Perhaps more de tailed free ion and J-J mixing calculations 
are needed to simulate upper levels satisfactorily. 
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CHAPTER -V 


SUMMARY 

In this thesis we have attempted an in-depth spectro- 
scopic study of Eu 3+ doped in LiYF 4 single crystal. When we 
undertook this investigation, no published work on the 
analysis of the spectroscopic data existed for this system. 
Subsequently Gorller et al [l] have reported an analysis of 
the fluorescence originating from the 5 D q and 5 D X multiplets 
of Eu 3+ . Our analysis of the LiYF^iEu 3 " 1 ’ data represents a 
somewhat different approach as compared to the one used by 
the authors of reference [.l]. Further we have considered 

several additional multiplets, not analysed so far. Gorller 

7 7 7 7 7 7 

et al have analysed F q , F^, F 2 , F^, F 4 and F^ multiplets 

5 5 ■=, 5 

whereas our analysis extends to D 0 , D^, D 2 and mul- 

plets as well. Moreover, unlike their work, we have re- 
corded polarized spectra and the polarization information 
was extremely useful in the symmetry assignment of levels. 

The absorption spectrum was obtained by using a IOOOW 

tungsten-halogen lamp and the fluorescence was excited by the 

457.9 nm and 514.5 nm lines of the Ar + laser and by a R-6G 

ring dye laser. The 457.9 nm line of Ar + laser non-reso- 

nantly excites the 5 D 2 multiplet of Eu 3+ and produces 

fluorescence from the ^D 2 , and 5 D Q multiplets. The 

5 5 

514.5 nm line produces fluorescence from the D x and D q 
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c* 

multiple ts. The D Q level can be excited selectively by 
tuning the R 6G dye laser to the 7 F^ ( ( ^ 1 ) transition. 

It can also be excited by a phonon assisted process when 

the dye laser is tuned to 17630 cm The excitation spectrum 

C »7 

was recorded by monitorinq the 3 D Q ( q)- F^ ( f^) transition at 
16840 cm -1 and scanning the dye laser from 16900 to 1780Q cm -1 .Tb 
excitation spectrum confirms the position of the D Q level 
at 17270 cm - '*' and has been useful in identifying Raman and 
vibronic lines in the fluorescence spectra. We have not 
attempted a serious analysis of the excitation spectrum 
because of the limited precision of these measurements 
(expected error in locating peaks is about 5-6 cm -1 ). The 
complete structure of the 5 D 1 and 5 D 2 mu 1 tip lets and partial 
structure of the 5 L 6 multiplet could be established from the 
absorption data. The 7 F o ~* 5 ° 0 and 5 °3 absorption could not 
be observed due to the forbidden nature of these transitions. 
5 D 2 - 7 Fj(J=0 to 6), 5 D 1 - 7 Fj(J=0 to 5) and 5 D Q - 7 Fj(J=0 to 
5)f luorescence has been recorded in it, cr and axial modes# 

All spectra were recorded at liquid nitrogen temperature# 

The Stark structures of the 7 F lf F 2 , F 3 and F^ have been 
established completely and the structure of the 7 F 5 could be 
established only partially. The analysis of the °2 “ F 6 
fluorescence has not been attempted because the lines belong- j 
ing to this group are rather closely spaced, weak and broad. 

For better understanding of \ - 7 F 6 fluorescence, it would j 
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be necessary to use a sample with hiqher Eu^ + concentration. 

Our experimentally inferred positions of the Stark levels 

7 . . 5 

belonging to the Fj(J=0 to 4) and D q ^ multiplets are in 
agreement (within +3 cm ) with those reported by Gorller 
et al. However, some disagreement exists for the structure 
of the multiple t. For example, our analysis of the 

experimental data would place the^ anci levels of the 
7 F 5 multiple t at 3786 cm" 1 and 3799 cm" 1 as compared to 
3795 cm -1 and 3807 cm” 1 reported by these authors. They 
have not observed the F^O^) levels at 3844 cm and 
4008 cm -1 which we have been able to establish. 

The CG's of the Stark components and hence the positions 

77777 7 

of the free ion levels of the F Q , F^, Fg» F^, F^, F^, 

-1 -1 

, ^D. and have been found to lie at 0 cm , 365 cm , 

1034 cm" 1 , 1921 cm” 1 , 2899 cm” 1 , 3948 cm , 17271 cm , 

19028 cm" 1 and 21502 cm" 1 respectively. The free ion 
Hamiltonian consisting of the spin-orbit, the residual 
Coulomb and the configuration interactions was solved numeri- 
cally by treating the radial quantities as adiustable para- 
meters in an effort to brina the calculated free ion energies 
of the nine multiplets in agreement with their experimentally 
inferred values. A systematic least squares analysis could 
not be carried out due to the large size of the free ion 
matrices. We had to resort to an approximate approach to 
perform this task. Only one parameter was varied at a time 
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till the minimum i n t ^ e r<rTUS< deviation was obtained for 
that parameter. This procedure was repeated for the remain- 
ing parameters. In this manner we could obtain an r.m.s. 
deviation of 24 crq 1 between the calculated free ion levels 

and the experimental CG ' s . 

The experimental crystal field splittings were obtained 
by taking differences of the observed Stark level positions 
with respect to the lowest Stark level of a given multiplet. 
Thus 23 such energy differences were obtained; 19 belonging 
to the Fj(J=0 to 5 ) a nd 4 to the ^Dj(J=0 to 2) multiplets. 

The multiplet was excluded from this analysis due to 

insufficient experimental data. 

The crystal field Hamiltonian was treated as a perturba- 
tion over the free ion Hamiltonian taking intermediate 
coupling wavef unctions as zeroth-order wavefunctions . Initi- 
ally crystal field calculations were performed neglecting 
J-J mixing effects i.e. crystal field matrix elements with 
J 4 J* were neglected. The least squares analysis was carried 
out for 19 energy differences and then for 23 energy dif- 
ferences, which gave r.m.s. deviations of 17.5 and 20.1 cm”'*' 

respectively. These deviations are considerably larger than 
— 1 

the 3 cm r.m.s. deviation reported by Gorller et al, who 

7 

considered J-J mixing effects within Fj manifold. If, 
however, we use crystal field parameters reported by these 
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authors, the r.m.s. deviations for 19 and 23 energy dif- 
ferences go up to about 23.2 and 24.8 cm”^ respectively. 

It was felt that large r.m.s. deviation may be due to our 
neglecting of J-J mixing effects. To estimate the contribu- 
tion of J-J mixing, we have carried out a limited J-J mixing 

calculation. Exact J-J mixing calculations are difficult to 

3+ 

perform for LiYF^jEu , because this needs large size matrices 
to be diagonalized. We have considered only two lowest 
multiple ts for each J (from J=0 to 6), and thus 6 energy 
matrices (corresponding to two multiplets and three irredu- 
cible representations) were constructed including matrix 
elements for J 4 J ' as well. These calculations gave very 

7 

good agreement for Fj multiple t but for upper levels agree- 
ment was poor. The least squares fit to the 19 energy dif- 
ferences gave a low r.m.s. deviation of 2.9 cm”'*' while that 

c 1 

for 23 levels (including the 3 Dj multiplets) was 13.1 cm - . 

7 

In fact the neglect of higher multiplets in case of the Fj 
manifold may be justified due to relatively pure nature of 

*7 

the 'F wavef unctions (96/ to 99/), but the same is not true 
for the higher levels. Other authors have encountered 
similar situation for Eu^ + in other hosts. The authors of 
ref [ 2 ] suggested the mixing of P states as a possible 
cause for the poor agreement for the Dj multiplets. This 
conjecture does not seem to hold, as we have used intermediate 
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3 

coupling wavefunc tions where contribution of P states are 
included. Though we do not have any satisfactory explana- 
tion for this poor agreement, we nevertheless feel that a 
more exact J-J mixing calculation involving other nearby 
multiplets should improve the situation but computational 
complexities have prevented us from performing such calcula- 
tions. It should be pointed out that we accounted for the 
free ion contributions in the J-J mixing calculations by 
adding to the diagonal matrix elements, the calculated free 
ion positions but in principle one should diagonalize the 
free ion and crystal field energy matrices simultaneously. 

We find that the order of calculated energy levels is 
in agreement with the experimental assignment based on 
polarization information in both the calculations (i.e. with 
and without J-J mixing effects). It may be mentioned that 
we always got the same values of the crystal field parameters 
irrespective of whatever values were used to initialize the 
minimization subroutine . 

We may conclude by saying that the LiYF 4 :Eu 3+ is a 
rather complex system and our effort to interpret the some- 
what incomplete data can be considered only a preliminary 
effort in that direction. Our description of the system is 
an improvement over the only other effort made so far. More 
detailed experimental and theoretical analysis especially 
for upper multiplets is needed for a better understanding 
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of this system. The role of phonons and other lattice 
effects vis-a-vis the observed temperature dependence of 
the positions of the spectral lines and the existence of 
efficient phonon assisted excitation effects are some of 
the other areas for future investigations. 


-xx- 
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APPENDIX-I 
TABLE AI-1 


LS TERMS OF 4f 6 CONFIGURATION 


S.No. 

S 

L 

(Wj Wp W3) 

(u x u 2 ) 

Nielson 
and Koster 
Notation 
(2S + l) Lr 

Senior! 
Quantum 
Number (’ 

1 . 

3 

3 

(100) 

(10) 

7F 

6 

2. 

2 

0 

(111) 

(00) 

5S 

4 

3 . 

2 

1 

(210) 

(11) 

5P 

6 

4 . 

2 

2 

■ (ill) 

(20) 

5D1 

4 

5 . 

2 

2 

(210) 

(20) 

5D2 

6 

6. 

2 

2 

(210) 

(21) 

5D3 

6 

7 . 

2 

3 

(111) 

( 10 ) 

5F1 

4 

8. 

2 

3 

(210) 

(21) 

5F2 

6 

9 . 

2 

4 

(111) 

(20) 

5G1 

4 

10. 

2 

4 

(210) 

(20) 

5G2 

6 

11. 

2 

4 

(210) 

(21) 

5G3 

6 

12. 

2 

5 

(210) 

(11) 

5H1 

6 

13. 

2 

5 

(210) 

(21) 

5H2 

6 

14. 

2 

6 

(111) 

(20) 

511 

4 

15. 

2 

6 

(210) 

(20) 

512 

6 

16. 

2 

7 

(210) 

(21) 

5K 

6 

17. 

2 

8 

(210) 

(21) 

5L 

6 

18 . 

1 

1 

( 110 ) 

(11) 

3P1 

2 

19. 

1 

1 

(211) 

(11) 

3P2 

4 
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20. 

1 

1 

(211) 

(30) 

3P3 

4 

21. 

1 

1 

(221) 

(11) 

3P4 

6 

22. 

1 

1 

(221) 

(30) 

3P5 

6 

23. 

1 

1 

(221) 

(31) 

3P6 

6 

24. 

1 

2 

(211) 

(20) 

3D1 

4 

25. 

1 

2 

(211) 

(21) 

3D2 

4 

26. 

1 

2 

(221) 

(20) 

3D3 

6 

27. 

1 

2 

(221) 

(21) 

3D4 

6 

28. 

1 

2 

(221) 

(31) 

3D5 

6 

29. 

1 

3 

(no) 

(10) 

3F1 

2 

30. 

1 

3 

(211) 

(10) 

3F2 

4 

31. 

1 

3 

(211) 

(21) 

3F3 

4 

32. 

1 

3 

(211) 

(30) 

3F4 

4 

33. 

1 

3 

(221) 

(10) 

3F5 

6 

34. 

1 

3 

(221) 

(21) 

3F6 

6 

35. 

1 

3 

(221) 

(30) 

3F7 

6 

36. 

1 

3 

(221) 

(31) A 

3F8 

6 

37. 

1 

3 

(221) 

(31) B 

3F9 

6 

38. 

1 

4 

(211) 

(20) 

3G1 

4 

39. 

1 

4 

(211) 

(21) 

3G2 

4 

40. 

1 

4 

(211) 

(30) 

3G3 

4 

41. 

1 

4 

(221) 

(20) 

3G4 

6 

42. 

1 

4 

(221) 

(21) 

3G5 

6 

43. 

1 

4 

(221) 

(30) 

3G6 

6 

44. 

1 

4 

(221) 

(31) 

3G7 

6 
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45. 

1 

5 

(110) 

(11) 

3H1 

2 

46. 

1 

5 

(211) 

(11) 

3H2 

4 

47. 

1 

5 

(211) 

(21) 

3H3 

4 

48. 

1 

5 

(211) 

(30) 

3H4 

4 

49. 

1 

5 

(221) 

(11) 

3H5 

6 

50. 

1 

5 

(221) 

(21) 

3H6 

6 

51. 

1 

5 

(221) 

(30) 

3H7 

6 

52. 

1 

5 

(221) 

(31) A 

3H8 

6 

53. 

1 

5 

(221) 

(31) B 

3H9 

6 

54. 

1 

6 

(211) 

(20) 

311 

4 

55. 

1 

6 

(211) 

(30) 

312 

4 

56. 

1 

6 

(221) 

(20) 

313 

6 

57. 

1 

6 

(221) 

(30) 

314 

6 

58. 

1 

6 

(221) 

(31) A 

315 

6 

59. 

1 

6 

(221) 

(31) B 

316 

6 

60. 

1 

7 

(211) 

(21) 

3K1 

4 

61. 

1 

7 

(211) 

(30) 

3K2 

4 

62. 

1 

7 

(221) 

(21) 

3K3 

6 

63. 

1 

7 

(221) 

(30) 

3K4 

6 

64. 

1 

7 

(221) 

(31) A 

3K5 

6 

65. 

1 

7 

(221) 

(31) B 

3K6 

6 

66. 

1 

8 

(211) 

(21) 

3L1 

4 

67. 

1 

8 

(221) 

(21) 

3L2 

6 

68. 

1 

8 

(221) 

(31) 

3L3 

6 

69. 

1 

9 

(211) 

(30) 

3M1 

4 
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70. 

1 

9 

(221) 

71. 

1 

9 

(221) 

72. 

1 

10 

(221) 

73. 

1 

11 

(221) 

74. 

0 

0 

(000) 

75. 

0 

0 

(220) 

76. 

0 

0 

(222) 

77. 

0 

0 

(222) 

78. 

0 

1 

(222) 

79. 

0 

2 

(200) 

80. 

0 

2 

(220) 

81. 

0 

2 

(220) 

82. 

0 

2 

(220) 

83. 

0 

2 

(222) 

84. 

0 

2 

(222) 

85. 

0 

3 

(220) 

86. 

0 

3 

(222) 

87. 

0 

3 

(222) 

88. 

0 

3 

(222) 

89. 

0 

4 

(200) 

90. 

0 

4 

(220) 

91. 

0 

4 

(220) 

92. 

0 

4 

(220) 

93. 

0 

4 

(222) 

94. 

0 

4 

(222) 


(30) 3M2 6 

(31) 3M3 6 

(31) 3N 6 

( 31 ) 30 6 

(00) 1S1 0 

(22) 1S2 4 

(00) 1S3 6 

(40) 1S4 6 

(30) IP 6 

(20) 1D1 2 

(20) 1D2 4 

(21) 1D3 4 

(22) 1D4 4 

(20) 1D5 6 

(40) 1D6 6 

(21) 1F1 4 

(10) 1F2 6 

(30) 1F3 6 

(40) 1F4 6 

(20) 1G1 2 

(20) 1G2 4 

(21) 1G3 4 

(22) 1G4 4 

(20) 1G5 6 

(30) 1G6 


6 
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95. 

0 

4 

(222) 

96. 

0 

4 

(222) 

97. 

0 

5 

(220) 

98. 

0 

5 

(220) 

99. 

0 

5 

(222) 

100. 

0 

5 

(222) 

101. 

0 

6 

(200) 

102. 

0 

6 

(220) 

103. 

0 

6 

(220) 

104. 

0 

6 

(222) 

105. 

0 

6 

(222) 

106. 

0 

6 

(222) 

107. 

0 

6 

(222) 

108. 

0 

7 

(220) 

109. 

0 

7 

(222) 

110. 

0 

7 

(222) 

111. 

0 

8 

(220) 

112. 

0 

8 

(220) 

113. 

0 

8 

(222) 

114. 

0 

8 

(222) 

115. 

0 

9 

(222) 

116. 

0 

9 

(222) 

117. 

0 

10 

(220) 

118. 

0 

10 

(222) 

119. 

0 

12 

(222) 


(40) A 1G7 6 

(40) B 1G8 6 

(21) 1H1 4 

(22) 1H2 4 

(30) 1H3 6 

(40) 1H4 6 

(20) 111 2 

(20) 112 4 

(22) 113 4 

(20) 114 6 

(30) 115 6 

(40) A 116 6 

(40) B 117 6 

(21) 1K1 4 

(30) IK 2 6 

(40) 1K3 6 

(21) 1L1 4 

(22) 1L2 4 

(40) A 1L3 6 

(40) B 1L4 6 

(30) 1M1 6 

(40) 1M2 6 

(22) INI 4 

(40) IN 2 6 

(40) IQ 


6 


APPENDIX-II 

DOUBLE MONOCHROMATOR GEM 1000 


The GDM-1000 double monochromator has two autocollimated 
monochromators, each with a concave mirror and a grating 
(651 lines/mm) in Littrow arrangement, placed one behind 
the other so that the dispersion of the gratings are added 
and coma aberration of mirrors annulled. This gives a sharp 
image of the entrance slit in the plane of the outlet slit 
and is free from coma. Other image aberrations are negligible 
due to the chosen focal length of the mirrors (f=1100 mm, 
aperture ratio 1:10.4) and the small angle of incidence of 
rays on the concave mirrors. A field flattening lens in 
front of the outlet slit is used to compensate for the 
curvature of the field arising from the concave mirrors and 
wavenumber dependence of the line curvature. The monochro- 
mator incorporates a built-in 25 Hz chopper and an 

automarking generator. 
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TABLE AII-1 
CALIBRATION CHART 

Source used for calibration : Ne discharge tube 
Spectral region 23000-17000 cm - ''' (lind order) 


Standard Wave No. 

(cm -1 ) .... . 

Observed on GDM 1000 
(cm -1 ) 

22737.9 

22738 

22580.8 

22580 

22281.2 

22279 

22037.3 

22036 

21475.7 

21475 

21256.7 

21256 

21236.6 

21236 

21231.1 

21231 

21222.1 

21221 

21207.4 

21207 

20881.5 

20881 

20715.3 

20714 

20471.2 

20469 

20173.4 

20172 

19850.1 

19849 

19436.6 

19435 

18759.0 

18759 

18722.8 

18721 

18715.1 

18714 

18516.6 

18515 

17396.5 

17395 

17347.8 

17347 

17086 . 8 

17084 

17001.3 

17000 



TABLE AII-2 


CALIBRATION CHART 

Source used for calibration : Ne discharge tube 
Spectral region 17350-11300 cm -1 (1st order) 


Standard Wave No. Observed on GDM 1000 

(cm" 1 ) (c m 1 ) 


17347.8 

17344 

17086.8 

17083 

17001.3 

16998 

16734.9 

16732 

16583.8 

16581 

16462.7 

16460 

16278.5 

16276 

16224.3 

16222 

16084.2 

16082 

15957.9 

15955 

15786.8 

15784 

15666 .6 

15664 

15619.5 

15617 

15153.9 

15152 

14973.9 

14972 

14431.1 

14429 

14219.9 

14218 

13939.4 

13938 

13802.3 

13801 

12047.7 

12046 

11936.6 

11935 


11767 

1 177 1.1 



11520 

11521.4 



11516 

11518.2 



11383 

11384.7 




TABLE All -3 
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P0 LAri 1 

ZA 'ION RESPONSE OF GDM-lOOO MONOCHROMATOR 

) 

Source u 

Spe c tral 

sed:1000*‘V Tungsten-! 

ronion: 25600 cm -11 

Halogen Lamp 

- 17500 cm -1 (Ilnd Order] 

Wave numb 
(cm"' 1 ' ) 

o r ( Arb 

H | | Slit 

(la) 

. uni ts ) 

E J. Slit 

(Im) 

I /I 
a rt 

25578 

1.4 

0.7 

2.0 

25178 

2.7 

1.5 

1.8 

24778 

4.6 

2.5 

1.84 

24378 

6.8 

4.7 

1.7 

23978 

8.7 

6.2 

1.4 

23578 

10.8 

8.6 

1.26 

23178 

1^.5 

11.0 

1.14 

22778 

13.7 

13.2 

1.04 

22578 

14.1 

14.1 

1.0 

22378 

14.2 

14.9 

0.95 

21978 

14.7 

16.4 

0.90 

21578 

14.8 

17.3 

0.86 

21178 

14.1 

17.5 

0.81 

20778 

13.2 

17.0 

0.78 

20378 

11.8 

16.0 

0.74 

19978 

10.1 

14.5 

0.70 

19578 

8.3 

12.6 

0.66 

19178 

6.6 

9.3 

0.71 

18778 

5.1 

8.2 

0.62 

18378 

3.9 

8.3 

0.47 

17978 

3.1 

7.9 

0.39 

17578 

2.2 

6 .6 

0.33 
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TABLE All -4 

POLARIZATION RESPONSE OF GDM-IOOO MONOCHROMATOR 


Source used: 1000 -W Tung s ten-Haloqen Lamp 
Spectral region: 17500-13000 cm - ' 1 ' (1st Order) 


Wavenumber 

( arb 

.units) 

I /I 

(cm ^ ) 

E | | Slit 

E 1 Slit 

a it 


(I ) 

(I ) 



a 

n 


17488 

17.3 

2.7 

6.41 

17288 

17.9 

2.8 

6.39 

17088 

18.3 

3.1 

5.90 

16888 

18.6 

3.5 

5.31 

16688 

18.8 

3.8 

4.95 

16488 

18.8 

.4.3 

4.37 

16288 

18.8 

5.0 

3.76 

16188 

18.6 

5.3 

3.51 

16088 

18.4 

6.4 

2.88 

15988 

17.8 

6.1 

2.92 

15888 

17.2 

5.5 

3.13 

15688 

15.8 

3.6 . 

4.39 

15488 

13.9 

1.9 

7.32 

15388 

13.0 

1.4 

9.29 

15288 

12.3 

1.2 

10.25 

15188 

11.4 

1.2 

9.5 

15088 

10.5 

1.2 

8.75 

14888 

8.8 

1.2 

7.33 

14688 

7.0 

1.2 

5 . 83 

14488 

5.5 

1.2 

4.58 

14288 

4.4 

1.2 

3.33 

14088 

2.7 

1.0 

2.70 

13888 

1.8 

0.7 

2.57 

13688* 

8.4 

5.4 

1.56 

13488* 

4.8 

3.3 

1 + 45 

13288* 

4.2 

2 .9 

1.45 

♦Recording 

parameters such 

as sensitivity and 

slit width were 


different . 



APPEND IX-III 


A. PROCEDURE FOR THE LEVEL ASSIGNMENT 

The assignment of transitions involving multiplets with 
J=0 (as one of the states) is straightforward and the posi- 
tions and symmetry (irreducible representations) of these 
states could be assigned easily using the polarization 
selection rules described in Chapters II and IV. All transi- 
tions in the absorption spectrum, the fluorescence transitions 
terminating on and those orginating from ^D o fall in this 

category. The level ass ionme n^. for fluorescence transitions 
originating from and terminating on multiplets with J ^0 
require more careful handling of the experimental data. The 
multiplets studied in this work are relatively well separated 
and in most cases it was possible to unambiguously associate 
a given group of fluorescence lines to transitions between a 
pair of multiplets. We will illustrate the procedure adopted 
in such cases by taking the example of ^-^3 fluorescence. 

The multiplets 5 D 1 and 7 F 3 have two (^+1^,4^ and flve 
( f^+2^+2(^ 4 ) components respectively. Thus a maximum of 10 
transitions are expected ignoring selection rules. There are 
9 observed lines from 17185 cm " 1 to 16984 cnT 1 which may 
belong to this group. Since we know the positions of both 
the components of \ from the analysis of other groups, we 
can construct the following table where the first column lists 
all the observed fluorescence lines and the second and third 
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columns give the difference of these from the positions of 

c 

Stark components. 


Observed transition 
wavenumber (cm -1 ) 

VQ.4> 

19021 cm -1 

Vr? , 

19043 cm -1 

17185 

1836 

1858 

17171 

1850 

1872 

17162 

1859 

1881 

17149 

1872 

1894 

17119 

1902 

1924 

17091 

1930 

1952 

17069 

1952 

1974 

17006 

2015 

2037 

16984 

2037 

2059 


From this table we notice that the lines at 17185 and 17162 
cm” 1 terminate on the same final level (1859+1 cm and 


their polarization behavior suggest them to be (T^ 4 ~ Hp and 
(fT-m transitions respectively (C.f. table 3.5 and 3.6) of 
Chapter III). Thus one of the levels of can be placed 

at 1859 cm” 1 . Similarly other and both levels of the 

7 F 3 can be placed at 2037 cm 1 , 1872 cm 1 and 1952 cm res- 
pectively. The remaining one transition at 17119 cm may be 
assigned to a level at 1902 or 1924 cm” 1 , but its polariza- 
tion behavior suggests it to be [ 3,4 transition. Since 

both 4 levels of 7 F 3 are already fixed, the only possi- 
bility left is to assign f[ at 1902 cm This assi 9 nment 

could be further confirmed by similar analysis of the F 3 

group. In this manner all the components of % could be 


identified 
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B. PROGRAMMING DETAILS 

1. FREE ION CALCULATIONS 

For the free ion calculations the energy matrices of 
the combined electrostatic, spin-orbit and configuration 
interactions are constructed. A typical matrix element is 

of the form 

E(I) = AE1(I) * E 1 + AE2 (i) * E 2 + AE3(l) * E 3 
+ ASO(I) *^ f + AAL(I) * a + ABE ( I ) * p 
+ AG A ( I ) * Y 1 

Where, the radial parts are denoted by their usual notations 
and the quantities multiplying them are the corresponding 
anqular pa rts * Since the energy matrix is real and symmetric, 
only one half of the matrix elements need to be calculated. 

For the Eu 3+ (4f 6 ) configuration, a total of 4228 matrix 
elements have to be calculated if the calculations are re- 
stricted to J values lying between 0 to 6 only. The angular 
parts of all the matrix elements were stored sequentially 
and were then broken into submatrices, one for each J. The 
order of the matrices corresponding to J = 0,1, 2, 3, 4, 5 and 6 
were 14,19,37,37,46,37 and 38 respectively. The angular 
parts of the Coulomb and configuration interactions were 
taken directly from the tables, whereas the angular parts of 
the spin-orbit interaction were calculated with the help of a 

Computations were done on DEC1090 system at Computer Centre, 
IIT Kanpur. 
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separate program using subroutines for calculating 3-J 
and 6-J symbols. The reduced matrix elements of were 

taken from tables. All these matrices were diagonalized 
using the 'EIGEN’ subroutine taken from the*IBM360 Scientific 
Subroutine Packaged The eigen values were printed out 

\ i- 

taking the energy eigen value of the lowest level as the 
zero of the energy scale . The eigen vector could be printed 
out whenever needed. 


2 * CRYSTAL FIELD CALCULATIONS 

The crystal field Hamiltonian in terms of the tensor 
operators is written as 


^crv^4-) — ^ " r u LLJ A vw 4 

cry 4 k=2, 4,6 0 0 k=4,6 4 4 


C k + E [B k (cj k) + cty) + iB* k 


'-4 

( C ^ k) - C ^ k) )] ...2 


A typical matrix elements of this Hamiltonian can be 
expressed as 


<f N SLJJ |E B k C^ |f N S ’ L’ J ’ J ’> = S ce , E F(k,q). ATOM(k) . 
z ' 44 z ss k,q 


... 3 


Where 

J-J 7 J k J’ 

F(k,q) = (-1) ( ) <f| |C Ck; | |f> 


. . .4 


“j, q J; 


ATOM(k) = (-D s+L,+J+k . vt (2J+D Tsn+m. 


fj J’ ki 

(L' l sj 


. <f N SL| |l/ k ^ 1 |f N SL’> 


• t • 5 


(k) 


The non-zero reduced matrix elements of U were taken from 
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the tables of Nielson and Koster were stored as such along 
with their SL values. The J,S and L values of the connecting 
states of the matrix elements were generated from the stored 
data within the computer program for all the matrices and 
were arranged in the sequential order compatible with the 
free ion wavef unction . The quantities F(k,q) and ATOM(k) 
were then calculated in that sequence and stored in a three 
dimensional array specifying the irreducible representation 
f£ j the (k,q) values and the sequence (Ni) of the matrix 
elements for a given J, as the three dimensions. Even zeros 
were included at this stage. 

The free ion wavefunctions were generated in a separate 

program by diagonalizing the free ion matrices and the 

products of the coefficients of the connecting | LS J> states 

(k) 

were combined with the corresponding U and summed over 
appropriately. They were subsequently multiplied by F(k,q)* s 
and then by the crystal field parameters to obtain the final 
crystal field matrix elements. 

2.1 CALCULATIONS WITHOUT J-J MIXING 

The crystal field energy matrices were constructed for 
each representation of a given J multiplet (i.e. the matrix 
elements between J and J' with J 4 J* were neglected). The 
resulting matrices were not more than 4x4 dimension but each 
matrix element was obtained summing over, all LS basis states 
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because of the intermediate coupling used for these calcula- 
tions. For example, the calculation of a single matrix 
element of a J wi th N | LS> basis states needs 7 ATOM(k) 's 
to be evaluated. The total number of were more than 

14x10 . The crystal field matrices were complex and hence 
were converted into real matrices of double the dimension- 
ality by a program written for this purpose and then 
diagonalized using the subroutine •EIGEN*. 

2 . 2 CALCULATIONS WITH J-J MIXIN5 

To include J-J mixing effects, the energy matrices were 
constructed for each representation for all "^Fj multiplets 
(J=0 to 6) and all 5d 0 ^ >2 3*^4 5 and mul ‘ t iP ie ' ts separate 

by now including the matrix elements with J ^ J* as well. 

This made calculations of the matrix elements more compli- 
cated and lengthy because now a given | JLS> basis of the 
wavefunction of a given J-multiplet were to be combined with 
all the |JLS> basis of the wavef unctions of other inter- 
fering multiplets. Though the resulting matrices did not 
exceed 13x13, the intermediate storage requirements increased 
immensely. For example, now the total number of ATOM(k) for 
the Fj multiplets alone was more than 10 . The program 
were modified not to store ATOM(k) and F(k,q) separately, 
but to be evaluated and used in the program itself. The 
calculated free ion eigen values were added to the diagonal 
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matrix elements of the crystal field energy matrices. 

These matrices were diagonalized in a similar fashion as 
for without J— J mixing. Several internal checks were made 
to ensure that no computational mistakes remain. 

3. minfun 

This program to minimize a function of several 
variables was written by W E Humphrey at Lawrence Radia- 
tion Laboratory [l]. This program was modified for some 
of its input requirements [ 2 ]. 

Theory 

This program uses the Ravine stepping method for 
finding the minimum of a function defined in the parameter 
space. A brief description of this program is given here. 
The variable names used here appear as a symbol in the 
MINFUN. 

Operation 1 ; The search procedure begins by computing the 
gradient GS of the function at the initial point 3^5 in 
the parameter space. The starting direction of search 
EXVEC is taken apposite to GB i.e. - GB. 

Operation 2 : A step is taken to point c such that 

XC = XB + (STEP .WT) EXVEC 6 

Where STEP is the step size and WT is the relative weight 
of the parameter. The value of the function FC and gradient 
GCf at this point are evaluated and a hyper-plane is defined 
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as to include GC and EXVEC. A unit vector PERP perpendicular 
to EXVEC and parallel to - GC + (GC. EXVEC) EXVEC is defined. 

Operation 3 : A step is taken along PERP from XC to another 
point XT = XC + (ALPHM.WT) PERP where ALPHM is determined 
in the programme and is of the order of twice the step 
size. Assuming the quadratic relation between function and 
ALPH, a point XA can be determined such that function FA is 
less than FC and FT where 

XA = XC + (ALPH.WT) PERP 7 

Operation 4 ; If function FA is truly a minimum (i.e. less 
than FC and FT) then the vectors EXVEC and XB are redefined 
as 

EXVEC = 8 

1 XA - XBj 

and 

XB = XB + (STEP .WT) EXVEC 9 

This completes a full step and control returns to 
operation 2 to begin the next step. If FA is not a new 
minimum, the control is returned to operation 3 with ALPHM 
half of its previous value. In the case that ALPH becomes 
more than twice the step size STEP, the previous parts of 
this step is skipped and the new step is taken along GC 
from XC. 



INPUT PARAME TERS 
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1. ISWS ( 


2. NPAR: 

3. NSTEP : 

4. STEP: 


5. EPST : 


6. X: 

7. WT: 

8. DIRIN: 


1-6): Each element of this array controls the 
corresponding sense switch and can have 
values 0 or 1. For details see Table 
(A III-l) 

Number of variable parameters n in the function. 

The maximum number of steps for the run. 

Magnitude of the step. If it is positive the 
search mode is used. If step is negative, 
convergence mode is used. If step is zero, it is 
set to +1.0 by the program. 

The tolerance. Normally it is zero. In the 
convergence mode if the function does not change 
by this much amount for each of the eleven 
consecutive steps, the run is terminated. 

Starting values of the parameters. 

The weight factors corresponding to the changes 
in the parameters. These must be non-zero. 

The initial direction for the first step. The 
individual elements of this vector may be of 
any size, as they are normalized by MINFUN. 
However, their signs are taken into account. 
These may be put as zeros, if no special direc- 
tion is required. 
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TABLE A-III-1 
SENSE SWITCHES 


Value of i 
option No. 

Comment 

1 

= 0 for maximum likelihood function 
= 1 for chi-squared function 

This affects only the error analysis. 

2 

not referenced-may be set to zero 

3 

Different minimizing procedures will 
be used by MINFUN, subject to certain 
tests, in order to give quicker 
convergence. 

4 

Error analysis is carried out at 
termination of run, including the 
error matrix and estimates of standard 
deviations of parameters. This option 
must not be used in the searching 
mode, as the error analysis is then 
carried out for the last point reached 
and not at the minimum. 

5 

Additional intermediate print-out is 
produced, describing each new maximum 
and minimum as it is discovered; and 
giving summary tables with graphs of 
the function and each parameter over 
each group of 50 steps. 

6 

= 0 for end of run 
= 1 if data follows 

Thus a blank card following all data 
ends the run. 
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9. FCNDATA : Data for the input to FCN subroutine, if any. 
3.1 MODIFICATIONS [2] 

Since no rigid criterion was available for the choice 
of WT S in help of MINFUN the program was modified so 
as to generate the WT'S internally and vector DIRIN was set 
to zero In the program. Thus no input data for WT and 
DIRIN was required by the modified version of the program. 
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